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Application  of  the  Esmigcoup  theory  of  operators  to  infirdte-di~ 
mensional  systems  is  an  exciting,  roeont  cleveloprcent  in  control 
theory.  In  tlie  Liid--19C0‘s  control  theorists  were  finding  scores  of 
n3\;  applications  for  the  rapidly  developing  ’'irodern  control  thecry"  of 
f inite“dii:;cnsiGri;l  systfiis.  It  is  f<air  to  say  that  tliis  new  theory 
was  met  with  eorisiderahle  rcsistanoe  and  skepticism;  and,  on  the  other 
hond,  Jt  \iC3  cer!  aiiPy  cv.rrr.ted  l:y  tour-  at  the  time.  Today,  two  dec¬ 
ades  l  et'.'",  applies. cd  o.:.>  of  i  odern  control  theory  are  still  being 
discover':-’.  Irof'cc].  its  uref  nl  r'ss  is  row  qiestior,'':d  by  ft-.',  cind 
tcchi  '.cr  :■  dovolop-'d  frer  f  1 1 j.t<>d:'.!i:on;.ional  redern  control  thcc'i'y  are 


rovi  .l  ly  nsyht  ir.  i:".:- v  <  rstt  .■  es  ru'orrd  the  v?crld,  .1  believe  scrii- 
gro  .'i;  li.i.  cry  to  br  c  cxtoi.sicn  of  t/iis  theory  to  inf  ir  i  ic-di- 

I  1  oys'i  1  an  corv;i.nccd  that  toe-  rudiiisnts  of 

thjf  tP  -■•y  vijJ  soor:  i  s>.ar:’a>-d  t  f-, ;•  ivo; '•  graduate  studeni  s  of 
can';'- 11,1  l,hec,'’y. 
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ANALYSIS  OF  THE  DYNAMIC  BEHAVIOR  OF  AN  INTENSE 


CHARGED  PARTICLE  BEAM  USING  THE  SEMIGROUP  APPROACH 


I.  Introduction 


Specific  Area  of  Research 

This  investigation  is  concerned  with  the  problem  of  controlling  a 
physical  system  which  is  most  naturally  described  by  a  set  of  partial 
differential  equations  (PDE).  The  many  successes  in  the  application 
of  the  "state  variable'  or  "modern  control  theory"  approach  to  systems 
of  linear  ordinary  differential  equations  (ODE)  have  led  many 
researchers  to  look  for  a  practical  extension  of  this  theory  to  accom¬ 
modate  systems  of  linear  PDE.  As  recently  as  1978,  however,  a  promi¬ 
nent  researcher  observed  (Russell,  1978:640):  "The  control  theory  of 
partial  differential  equations  has  followed  right  on  the  heels  of  that 
for  ordinary  differential  equations,  but  with  slower  and  heavier 
tread . 

Models  for  many  physical  systems  can  be  brought  into  the  form  of 
an  abstract  C a u c Ii v  nr ob len.  Let  X  be  a  Banach  space,  and  suppose  A 
is  a  linear  operator  from  a  subset  of  X  into  X  .  If  the  domain  of  A 
is  dense  in  X  ,  then  the  equation 

u  ("  t  )  =  A  u  (  t  )  t  '■  0 


I-l 


and  is  denoted  by  5f  .  Existence  of  the  Gateaux  derivative  of  F  at 

X 

X  does  not  imply  continuity  of  the  operator  F. 

On  the  other  hand,  one  can  generalize  the  derivative  of  an 
operator  in  a  manner  which  mimics  the  "differentiability  implies 
continuity"  property  of  the  usual  derivative.  Suppose  for  some 

A 

x-;X.  ,  an  open  subset  of  V(F)  ,  there  exists  a  5F^£B(X,Y)  such 

that 

A 

lim  |jF(x+h)  -  F(x)  -  6F^(h)||Y 

The  operator  5F^  is  termed  the  Frechet  derivative  of  F  at  x,  and 

the  existence  of  this  derivative  implies  continuity  of  F  at  x 

A 

(Luenberger,  1969;  173).  Furthermore,  existence  of  oF^  implies 

existence  of  oF  and  the  two  are  equal  in  this  case. 

X 

The  Gateaux  and  Frechet  derivatives  are  often  used  to  construct  a 
linear  approximation  of  a  nonlinear  operator.  The  procedure  is 
analagous  to  the  familiar  first-order  Taylor  series  linearization 
techniques  for  a  real  function  of  a  real  variable. 

Consider  next  a  function  u:I-^X,  where  I  is  an  interval 
(possibly  infinite)  of  the  real  line,  and  X  is  a  Banach  space.  If 
the  Frechet  derivative  of  u  at  t  ^  £  I  exists,  then  this  operator  is 
called  the  strong  derivative  of  u  at  t^.  For  this  special  case,  the 

A 

cumbersome  Frechet  derivative  notation  is  replaced  with  the 

d  .  “ 

usual  differentiation  symbols  t  )  or  u(t^)  . 
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D^f  = 


'  k  I 

3  '  '  f 


kl  ^  k; 


d  X  , 


^  k 
ox  ri 
n 


where  k=(k^  ,k,  ,  .  .  -k^)  ,  the  k^  being  nonnegative  integers,  and, 

n 

!k'  =  Zk 

i=l 

Unless  stated  otherwise,  use  of  the  usual  differentiation  symbols  D 

""k 

or  ^ —  indicates  a  generalized  derivative.  Various  functional 

0  X, 

k 

analysis  texts  cover  generalized  derivatives  in  detail  (also  known  as 
distributional,  and,  more  generally,  as  weak  derivatives).  See 
(Curtain  and  Pritchard,  1977:  136-138;  Yosida,  1968:  48-52),  for 

example. 

Consider  now  an  operator  F  (not  necessaily  linear)  with  domain 
P(F)  a  subspace  of  a  normed  linear  space,  X,  and  range  contained  in 
a  normed  linear  space  Y  .  Two  generalizations  of  the  derivative  of  a 
real  function  of  a  real  variable  are  possible,  where  an  appropriate 
topological  generalization  of  R  is  assigned. 

First,  consider 


1 i m  F ( x+h V )  -  F ( x ) 

h-^0  h 

where  x  ,  veP(F)  ,  and  hcR.  if  the  limit  exists  for  every  veP(F), 
then  the  operator  F  is  said  to  be  Gateaux  differentiable  at  x  .  In 
this  case,  the  limit  above  defines  a  unique  element  in  Y  for  every 
v-:D(F)  .  This  mapping  is  called  the  Gateaux  derivative  of  F  at  x 
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denoted  bv  C(X,Y)  ,  or  by  C(X)  if  X=Y  .  For  examples  and  further 
discussions  of  closed  operators  the  reader  is  referred  to  (Curtain  and 
Pritchard,  1977:  45;  Belleni-Morante,  1979:  60-63;  Taylor  and  Lay, 

1980:  208-217). 

Derivatives.  Various  generalizations  of  the  usual  derivative  of 
a  real  function  of  a  real  variable  exist,  depending  on  the  topological 
properties  assigned  to  the  underlying  spaces.  Three  specific  types  of 
derivatives  are  of  use  in  the  application  of  semigroup  theory: 
(1)  generalized  derivatives,  (2)  Gateaux  derivatives,  and  (3)  Frechet 
derivatives.  The  Frechet  derivative  of  a  function  whose  domain  is  an 
interval  of  the  real  lir^  is  known  as  a  strong  derivative.  Since  the 
strong  derivative  is  frequently  used  in  semigroup  theory,  it  is  also 
discussed  below. 

Let  (il)  denote  the  set  of  all  functions  'P  which  are 

continuous,  have  continuous  partial  derivatives  of  any  order,  and 

which  have  support  bounded  and  contained  in  Q,  an  open  subset  of  , 

The  generalized  derivative  of  any  function  feL^  (f^)  (  f  only 

1  o  c 

"locally"  belongs  to  L^  (P. )  ,  i.  e. ,  f  is  defined  on  0.  and  is  in 

(K)  for  every  Lebesgue  measurable  set  K  whose  closure  is 

contained  in  -  ) ,  if  it  exists,  is  defined  to  be  the  function  g  such 
that 


f(x)D*^J  (x)dx  =  (-1)*^  /  g(x)i  (x)dx 


,  f 

g  (x 


for  all  ? C  (f. )  .  The  differentiation  operator  D  is  defined  as 
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the  following  two  statements  hold  (Naylor  and  Sell,  1982;  240): 


(i)  If  T  is  continuous  at  any  point  x  eX  ,  then  it  is 
continuous  at  every  xeX, 

(ii)  T  is  bounded  if  and  only  if  it  is  continuous. 

Although  bounded  linear  operators  are  simpler  to  analyze, 
unbounded  linear  operators  frequently  appear  in  applications.  The 
"derivative"  operator,  for  example,  is  often  unbounded,  depending  upon 
the  domain  and  codomain  chosen  for  a  specific  model.  In  some  cases  an 
unbounded  linear  operator  enjoys  properties  similar  to  those  of  a 
continuous  one,  in  which  case  it  is  termed  a  closed  operator.  The 
definition  of  a  closed  operator  is  often  stated  in  terms  of  its  graph, 
but  an  equivalent  and  more  practicail  definition,  in  the  context  of 
metric  spaces,  is  the  following: 

Definition  2.2  (Closed  Operator) 

Let  T  :  D  ( T  )CX->-Y  be  an  operator  with  x,  Y  Banach  spaces. 

Suppose  is  a  sequence  in  P(T)  with  the  properties 

(1)  X  ^x 

n 

(ii) 

The  operator  T  is  closed  if  xeP(T)  and  Tx=y  for  every 

such  sequence  in  V(T)  , 

The  set  of  all  closed  lineair  operators  defined  on  a  subset  of  the 
Banach  space  X  and  with  range  contained  in  the  Banach  space  Y  is 
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X  are  c;:p.’ c  ■  ly  ui'-r'id  ir:  a-yl  ic ;  1  i  oi:;. .  Tho  J  iyyXc^lly 

ropvv,:' ^nls  vaUiar,  of  wluAe  l),e  E.riaol,  X  la  iv'oat  ona.a  a 

spac.o  cr  fiucJ.cns  (r.y.,  I/'  (  )  o-  il'(.)).  Suriaae  u  aaot.ea 

such  i.ipi  isy  (1.  n.  ,  r.  rICh'X  ),  then  in  aocnrdc.ncc  \.  ith  Iv.fl/  itiOM 


,  1,  u  is  c:iiit:Ltu-ci;.c  s-'o  L^f  T  if  f oi’  evti-y  u>0  uhci 


C'  o?'.:.  c  r 


ruch  lh;:t 


'>(•  Jllx 


jatl  ervinf 


it  ~  1,1  •  ^ 


hi  s.  :X',,v  enp  li  ' .  i  ,  .  lyninys  from  cn  ii/icrxil  I  of 


(;.  i:  ■; 


C'  U  •  •  J  C  1':  J.  Oj  I  i  I  0 


eX  r'(')  J 


.1  "  t'-.  ‘  ;  c  rrj.  O',  rir  i  Li  o;i  of 


;  rr  ■' 


J  :  ;  '  x.,  1  :  1  -  D  (X)  rt 


’  :  1.  Ic  ,  ;  ; 


-rid 


Linear  Operators  on  Banach  Spaces 


There  are  several  key 


concepts  involving  linear  operators  whose  domains  and  codomains  are 
subsets  of  normed  spaces,  A  bounded  linear  operator  T;X-+Y  ,  X,  Y 
Banach  spaces,  is  one  for  which  a  nonnegative  real  number  K  exists 
such  that 


Klifll 


X 


for  all  f£X.  The  set  of  all  bounded  linear  operators  from  X  into  Y 
is  itself  a  Banach  space  and  is  denoted  by  B(X,Y)  ,  or,  if  X=Y  ,  by 
B(X)  .  The  Inflmum  of  the  set  of  all  constants  K  satisfying  the 
above  inequality  is  the  norm  of  T  on  the  Banach  space  B(X,Y)  or 
B(X)  . 


Continuity.  The  notion  of  continuity  of  a  function  is 
fundamentauL  in  semigroup  theory  analysis.  The  following  definition  is 
sufficiently  general  for  subsequent  discussions: 


Definition  ZA.  (Continuity) 

Let  X  ,  Y  be  normed  linear  spaces,  and  suppose  F 
represents  a  function  from  a  subset  D  of  X  into  Y  — 
i.  e. ,  F:DCX^Y  .  F  is  said  to  be  continuous  at  the  point 
in  D  if  for  every  real  number  e>0  there  exists  a  real 
number  5  such  that 


i|F(x)  -  F(x„)I|y  "  " 
for  all  xeD  satisfying 


X  - 


^olx  ^ 


Mappings  from  an  interval  I  on  the  real  line  into  a  Banach  space 
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Various  linear  spaces  are  used  in  this  work.  The  set  of  real 

numbers  and  the  set  of  complex  numbers  are  symbolized  by  R  and  C  , 

respectively,  and  these  are  the  only  scalar  fields  used.  The  symbols 

R^  and  denote  n-fold  Cartesian  products  of  the  linear  spaces  R 

and  C  (with  the  usual  addition  and  scalar  multiplication 

definitions).  The  letters  I  and  are  used  to  mean  an  interval  of 

R  or  R^  ,  respectively,  either  finite  or  infinite  —  i.  e. , 

I  =  (a,b)CR  ,  and  2=  {  x£  :  x=  (x  ,  •  •  •  x  )  ,  a  .  <  x  .  <  b  .  ,  i=  1 .  .  .  nlC 

with  a  ,  a^eR  or  and  b,b^eR  or  +<» .  Occasionally  Cartesian 

products  of  linear  spaces  are  denoted  by  the  product  symbol,  H  . 

Specifically,  letting  ^  linear  spaces,  the 

n 

Cartesian  product  of  these  spaces  is  written  as  11  The  most 

common  function  spaces  used  in  this  report  are  the  Lebesgue  and 
Sobolev  spaces,  L^(r2)  and  .  The  spaces  consist  of 

(equivalence  classes  of)  functions  f  such  that  I  is  integrable  in 
the  Lebesgue  sense.  The  Sobolev  spaces,  H'^(f^),  consist  of  the  sets 
of  functions  f  whose  generalized  derivatives  (discussed  below)  up  to 
eind  including  order  q  are  in  L^(f2)  .  For  further  discussion  of 
Lebesgue  and  Sobolev  spaces,  the  reader  is  referred  to  (Hoyden,  1968: 
Ch  6)  and  (Yosida,  1968:  55),  respectively.  The  Sobolev  spaces  are 
Hilbert  spaces  for  all  integer  q^O ,  as  is  L^(n)  ,  and  L^(J2)  is  a 
Banach  space  for  all  integer  p^l.  The  norm  of  a  function  f  in  these 
normed  spaces,  or  any  other  normed  space,  is  symbolized  by  ||f||j^ 
where  X  represents  the  space,  or  by  ]|^|1’  ^  is  clear  which  space 


is  intended. 


II.  Pertinent  Results  from  Operator  Semigroup  Theory 


Introduction 

Some  known  results  from  the  semigroup  theory  of  operators 
(hereafter  referred  to  as  semigroup  theory)  are  now  presented.  The 
theory  has  been  rigorously  developed  by  Hille  and  Phillips  (1957), 
More  recent  texts  have  been  written  which  emphasize  the  practical 
aspects  of  the  theory  (Butzer  and  Berens,  1967;  Bell eni-Mor ante,  1979; 
Curtain  and  Pritchard,  1977;  Curtain  and  Pritchard,  1978;  Davies, 
1980;  Fattorini,  1983;  Pazy,  1983;  Walker,  1980),  The  intent  of  this 
chapter  is  to  state  notation,  definitions,  and  specific  results 
pertaining  to  semigroup  theory  and  the  abstract  Cauchy  problem  which 
are  relevant  to  an  analysis  of  the  beam  models  developed  in  the 
following  chapter. 


Fundamental  Notation  Definitions 

Functions  and  Spaces.  Let  A,  B  be  arbitrary  sets.  The  notation 
f:A->B  is  used  to  denote  a  function  f  with  dom£d.n0(f)  equail  to  A  , 
and  range  R(f)  a  subset  of  the  codomain,  B.  The  term  operator  is 
used  to  denote  any  function  whose  domain  or  codomain  (or  both)  is  a 


apace  of  functions 


these  readers  should  follow  the  remaining  chapters  with  little  diffi¬ 
culty. 

Chapter  III  introduces  various  dynamic  models  of  charged  particle 
beams.  The  sophisticated  "microscopic"  models  are  presented  first, 
and  a  linearization  is  performed  in  order  to  bring  this  class  of  mod¬ 
els  into  the  abstract  Cauchy  problem  form.  "Macroscopic"  descriptions 
are  then  discussed  in  general,  and  a  linear,  single  degree  of  freedom 
model  is  derived.  Finally,  a  tractable  model  is  developed  in  detail 
in  order  to  illustrate  semigroup  theory  techniques  analytically. 

An  analytical  solution  of  the  "electrostatic  approximation  model" 
is  thoroughly  developed  in  Chapter  IV.  Various  simplifying  assump¬ 
tions  are  introduced  in  Chapter  III  in  the  development  of  this  model 
which  would  not  necessarily  be  required  if  a  numerical  solution  were 
sought.  It  is  considered  far  more  useful  from  a  researcher's  point  of 
view,  however,  to  develop  a  closed-form  solution  to  the  electrostatic 
approximation  model  thoroughly  than  to  resort  to  a  numerical  solution 
of  a  more  complicated  model. 

A  summary  of  dissertation  research  results  is  presented  in  the 
concluding  chapter,  along  with  some  suggested  further  areas  of 


research 


nuclear  fusion  research  has  been  treated  in  a  manner  similar  to  that 


herein  (cf.  Wang  and  Janos,  1970).  The  plasma  confinement  problem 
differs  considerably  from  the  beam  dynamics  problem,  however.  The 
plasma  in  Wang's  work  is  assumed  to  be  neutral,  while  a  charged  parti¬ 
cle  beam  is  a  nonneutral  plasma.  Furthermore,  the  configuration  of 
the  plasma  confinement  problem  does  not  at  all  match  that  of  the  beam 
problem,  where  a  large  velocity  field  in  one  direction  is  assumed. 
Nonetheless,  the  starting  point  for  both  plasma  confinement  in  Wang's 
paper  and  the  beam  dynamics  analysis  in  this  dissertation  is  the 
Vlasov-Maxwell  system  of  equations. 

QYsrview 


Three  major  topics  are  presented  in  the  sequel:  (1)  a  summary  of 
relevant  mathematical  concepts,  (2)  a  description  of  various  mathe¬ 
matical  models  of  the  dynamics  of  a  charged  particle  beam,  and  (3)  an 
illustration  of  the  theory. 

The  purpose  of  Chapter  II  is  twofold.  First,  it  provides  readers 
with  functional  analysis  and  operator  semigroup  theory  in  their  back¬ 
grounds  a  summary  of  notation,  definitions,  and  results  in  these 
areas.  Second,  readers  of  this  chapter  with  finite-dimensional  modern 
control  theory  in  their  backgrounds  are  provided  a  glimpse  of  how  the 
finite-dimensional  theory  generalizes  to  the  infinite-dimensional  the¬ 
ory.  For  example,  a  real  matrix  operator  of  order  a  is  discussed  as  a 
special  case  in  the  subsection  "Some  Familiar  Operators."  Armed  with 
these  insights,  and  intuition  provided  by  finite-dimensional  theory. 
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primary  goal  of  this  dissertation  is  to  advance  the  development  of 
semigroup  theory  techniques  by  attacking  a  specific  initial  value 
problem:  the  dynamic  behavior  of  an  intense  charged  particle  beam. 

Intense  beams  of  charged  particles  are  beginning  to  be  used  in  a 
wide  variety  of  applications  (Septier,  1983:  xii).  The  dynamic  behav¬ 
ior  of  such  beams  is  quite  complex  because  electromagnetic  fields  are 
affected  by  not  only  the  positions  of  the  particles,  but  by  their 
velocities  as  well.  Frequently  the  Vlasov-Maxwell  system  of  PDE  is 
chosen  as  a  starting  point  for  analysis  of  a  collection  of  charged 
particles.  Simplifying  assumptions  are  often  appropriate,  but  the 
resulting  models  are  generally  systems  of  PDE  also.  Analysis  of  the 
dynamic  behavior  of  intense  charged  particle  beams  is  an  excellent 
choice,  then,  for  an  application  of  semigroup  theory  since  (1)  such 
beams  are  useful,  and,  (2)  models  of  these  beams  are  inherently  dis¬ 
tributed  parameter  systems  of  equations. 

This  dissertation  establishes  a  framework  for  analyzing  the  beam 
dynamics  problem.  In  the  figure  on  page  1-5  the  basic  problem  is 
divided  into  two  sub-problems:  (1)  the  control  problem,  which  is  con¬ 
cerned  with  modifying  the  dynamic  behavior  to  achieve  some  desired 
state,  and  (2)  the  observation  problem,  which  is  concerned  with  deter¬ 
mining  the  present  state  of  the  beam.  The  foundation  laid  in  this 
work  is  original  and  should  serve  to  direct  and  organize  beam  dyamics 
research  in  the  future. 

Some  articles  exist  in  the  literature  which  are  related  to  this 
research.  For  example,  the  plasma  confinement  problem  associated  with 


A 


+  8  (x ) 


m 

k=  1  k 

and  where  A  (x)  =  fa^.(x)}  and  8(x)  =  f8..(x)}are  n*"^ -order 

1 J  1  j 

matrix  functions  defined  in  R"’. 

It  has  been  shown  that  ordinary,  partial,  stochastic,  and  delay 
differential  equations  can  all  be  accommodated  by  the  application  of 
semigroup  theory  to  initial  value  problems  on  a  Banach  space  (Curtain 
and  Pritchard,  1978:  Ch  8).  Bell eni-Mor ante  (Belleni-Morante,  1979: 
Ch  8-13)  discusses  in  detail  the  following  specific  problems:  heat 
conduction  in  rigid  bodies,  one-speed  neutron  transport,  kinetic 
theory  of  vehicular  traffic,  the  telegraphic  and  wave  equations,  the 
one-dimensional  Schrodinger  equation,  and  stochastic  population 
theory.  Additionally,  Markov  processes  were  studied  from  the 
semigroup  theory  point  of  view  by  Hille,  Yosida  and  Feller  in  the 
early  1950*s  (Fattorini,  1983:  98).  These  examples,  and  many  others 
that  can  be  found  in  the  recent  literature,  illustrate  the  wide 
variety  of  physical  problems  that  can  be  formulated  and  analyzed 
within  the  context  of  semigroup  theory. 

This  diversity  of  applications  is  encouraging,  but  far  more 
practical  applications  are  needed.  Fattorini  (Fattorini,  1983:  xx) 
states,  "Nowadays,  many  volumes  devoted  ...  to  the  treatment  of 
semigroup  theory  exist...  In  contrast,  accounts  of  the  applications  to 
particular  partial  differential  equations  ...  are  scarcer..."  This 
suggests  that  more  applications  should  be  attempted  in  order  for  the 
theory  to  develop  into  a  practical,  working  body  of  knowledge.  The 
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along  with  an  initial  condition, 


u(0)  =  u° 

is  termed  an  abstract  Cauchy  problem. 

Analysis  of  the  abstract  Cauchy  problem  can  be  performed  with  the 
semigroup  theory  of  operators.  This  approach  has  many  parallels  with 
the  modern  control  theory  approach  to  systems  of  linear,  time 
invariant,  first-order  ordinary  differential  equations: 

x(t)  =  Ax(t)  t^O 

where  A  is  an  n  by  n  real  matrix.  For  example,  the  state  transition 

A.  t 

matrix,  e  ,  for  such  a  system  of  ODE,  is  an  element  of  a  semigroup 
of  operators  generated  by  A  ,  where  t^O  .  Another  parallel 

exists  in  that  the  semigroup  theory  onphasizes  spectral  properties  of 
the  operator  A  in  the  abstract  Cauchy  problem.  This  is,  of  course, 
analogous  to  the  modern  control  theory  emphasis  on  the  eigenvalues  and 
eigenvectors  of  the  matrix  A  .  These  parallels  provide  a  compelling 
case  for  considering  an  appropriate  extension  of  modern  control  theory 
to  be  analysis  of  the  abstract  Cauchy  problem  through  the  semigroup 
theory  of  operators.  This  point  of  view  is  adopted  in  the  present 
work. 

Only  llneeir  systems  of  partial  differential  equations  are 
considered  herein.  In  fact,  all  models  are  of  the  form 

■^w(x,t)  =  Aw(x,t) 

where  w(x,t)£R'^  x=(x  ...x  )eR"*,  A  is  given  by 
^  1  m  '  . 
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Spectral  Analysis  Definitions.  It  is  well  known  that  the 


eigenvalues  and  eigenvectors  in  a  finite-dimensional  system  of  linear, 
first-order,  time- invariant,  differential  equations  are  instrumental 
in  an  analysis  of  such  a  system.  In  infinite-dimensional  systems  the 
eigenstructure  is  equally  important. 

Let  A  :V  (A)->-X  ,  V(A)CX  ,  be  a  linear  operator  with  X  a  Banach 
space.  The  set  of  all  complex  numbers  can  be  partitioned  into  two 
subsets  according  to  whether  XI-A  satisfies  the  following  three 
conditions  for  AeC  (Yosida,  I968:  209;  Curtain  and  Pritchard,  1977: 
163,  164;  Naylor  and  Sell,  1982:  414-429): 

(i)  (Al-A)~'  exists 

(ii)  (Al-A)“^  is  continuous 

(iii)  the  range  of  A  I- A  is  dense  in  X 

The  set  of  all  AeC  such  that  these  conditions  are  met  is  called  the 
resolvent  set,  while  the  set  of  ail  other  complex  numbers  is  called 
the  spectrum.  The  resolvent  set  and  spectrum  are  denoted  by  P(A) 
and  o  (A)  ,  respectively. 

The  spectrum  of  a  linear  operator  A  defined  on  a  subset  of  a 
finite-dimensional  space  E  ,  with  range  contained  in  E,  consists  of 
only  those  XeC  such  that  Al-A  is  not  injective.  Similarly, 
Al-A  can  fail  to  be  injective  for  some  values  of  A  ,  but,  unlike 
the  finite-dimensional  case,  the  spectrum  may  contain  other  complex 
numbers.  In  fact,  there  are  three  disjoint  subsets  of  a  (A)  ,  The 
point  spectrum  consists  of  those  Agc  for  which  Al-A  is  not 
injective.  The  continuous  spectrum  is  made  up  of  those  A  for  which 
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(XI-A)~  exists  but  is  not  continuous,  and,  for  which  the  range  of 
Xl-A  is  dense  in  X.  Finally,  the  residual  spectrum  consists  of 
those  X  for  which  (XI-A)~^  exists  and  is  continuous,  but  such 

that  the  range  of  XI-A  is  not  dense  in  X  , 

Let  the  notation  R(z,A)  denote  the  operator  (zI-A)  *  for 

any  zco  (A)  .  The  following  two  facts  are  established  in 

(Bell erJL-Mor ante,  1979:  62,63): 

(i)  If  A  is  a  closed  linear  operator  (AeCCX)  ),  A^  is  a 
bounded  linear  operator  (  A^eB(X)),  and  if  the  domain 
of  Aj^  contains  the  domain  of  A  ,  then  A+A^£C(X) 

(ii)  If  for  any  ZgSC  ,  R  (z  ^  ,  A)  e  B  (X)  ,  then  AeCCX)  . 


These  two  facts  are  used  frequently  in  practical  applications  of  the 
semigroup  theory  of  operators.  For  example,  see  (Belleni-Morante, 
197y:  179)  where  the  first  fact  is  used  in  proving  an  important 

perturbation  theorem. 

The  set  of  closed  linear  operators  is  frequently  partitioned  in  a 
manner  which  simplifies  semigsroup  theory  discussions.  The  four 
classes  of  interest  are  denoted  by  G(1,S),  G'(',S),  G(M,S)  ,  and 
G'(M,3)  and  are  defined  as  follows  (Belleni-Morante, 

1979:140,141 ,145): 

Definition  2.3  (  G -Classes) 

Let  AeC(X)  ,  D(A)  dense  in  X,  zeC,  and  ;  =  Re(z). 

Then  A  is  in  the  class 

(i)G(l,3)  if  {z:r>S}CD(A)  and  ||r(2,A)|<t^ 

for  all  z  such  that  ?;>3 

(ii)  G'(l,3)  if  {z;|;|>3}Co(A)  and  ||r  (  z  ,  A  )||  £  i  ^ 


for  all  z  such  that  li;|>S 

(ill)  G(M,S)  if  {z:C>S}Cp(A)  and  for  any 
integer  j  =  1 , 2  ,  .  .  . 


for  all  z  such  that  ^>B 

(iv)G'(M,S)  if  {  z  ;  I  C  I  >B}Cp  (A)  and  for  any 
integer  j  =  1 , 2  ,  .  . . 


R(z  ,A)^ 


M 


(U1  -8) 


for  all  z  such  that  |c!>6 


The  various  mathematical  symbols  which  have  appeared  in  this 
section  are  summarized  in  Appendix  A.  With  these  fundamental 
definitions  and  results  in  mind,  attention  is  now  turned  to  the 
abstract  Cauchy  problem. 


The  Abstract  Cauchv  Problem 

Mathematical  models  are  frequently  developed  to  predict  the 
dynamic  behavior  of  certain  variables  in  a  physical  system.  In  many 
cases  the  model  is  finite-dimensional  and  one  is  interested  in  knowing 
what  values  in  R  each  vau'iable  assumes  at  any  given  time.  In 
distributed  systems,  however,  the  variables  of  interest  can  be 


elements  of  a  function  space  at  each  instant  of  time.  Quite  often  a 


2.6  (Well  Posed) 


The  abstract  Cauchy  problem  (2.1),  (2.2)  is  well  posed  in 
if  the  following  two  conditions  are  satisfied: 

(i)  Existence  of  solutions  for  sufficiently  many 
initial  data:  There  exists  a  dense  subspace  D 
of  X  such  that,  for  anv  u°eD  ,  there  exists 
a  solution  of  the  abstract  Cauchy  problem. 

(ii)  Continuous  dependence  of  solutions  on  their 
initial  data:  There  exists  a  nondecreasing, 
nonnegative  function  C(t)  defined  in  t>0 
such  that 


u(t) 


<  C(t) 


u(0)| 


(2.3) 


for  any  solution  of  the  abstract  Cauchy  prob¬ 
lem. 


These  requirements  are  similar  to  those  generally  deemed  essen¬ 
tial  in  order  for  a  mathematical  model  to  correspond  to  physical  real¬ 
ly  ity  (e.g.,  see  Courant  and  Hilbert,  1962:  227):  (1)  existence  of  solu¬ 

tions,  (2)  uniqueness  of  solutions,  and  (3)  continuous  dependence  of 
the  solution  on  the  initial  data.  For  instance,  the  well  posed  Cauchy 
problem  has  the  existence-of- solution  property  for  a  particular  set  of 
initial  conditions.  (However,  a  solution  is  not  always  guaranteed  to 
exist  for  every  u^eX,  but  only  for  every  u in  a  dense  subset  of  X.) 
Furthermore,  equation  (2.3)  ensures  that  any  solution  of  a  well  posed 
abstract  Cauchy  problem  is  unique.  To  demonstrate  this,  let  v,  w  be 
solutions  of 


^  u(t)  =  Au(t) 


(t>0) 


n 


f> 


and  consider  the  vector  v-w 


Clearly,  v-w  is  also  a  solution 


and,  since  (v-w)  (0)=0  ,  we  have  from  equation  (2,3) 


( v-w) ( c ) 


C(t)  0 


0 


from  which  it  follows  that  v=w  .  Heurist ical ly,  the  third 

requirement  is  that  two  initial  conditions  which  are  "close"  to  each 
other  should  yield  solutions  which  are  also  "close."  The  second  po- 
sedness  condition  ensures  this  continuous  dependence  of  solutions  on 
initial  data. 

Any  solution  of  a  well  posed  abstract  Cauchy  problem  with  initial 
condition  lying  in  D  (the  set  referred  to  in  the  first  posedness  con¬ 
dition)  uniquely  defines  an  operator  S(t);D->-X  as 

u  ( t )  =  S (t ) u° 

for  t  ^  0,  with  u(0)  =  u®  .  Furthermore,  S  ( t )  is  necessarily  a 
linear,  bounded  operator  in  D  (by  the  linearity  of  A  in  equation 
(2.1),  and  by  the  second  posedness  condition,  respectively)  and,  as  D 
is  dense  in  X  ,S(t)can  be  extended  to  all  of  X  ,  The  operator-valued 
function  S  is  called  the  propa2ator  for  the  solution  of  the  well 
posed  abstract  Cauchy  problem. 

Well  posedness  of  the  abstract  Cauchy  problem  supports  a  notion 
of  "solution"  for  any  u®eX  .  Indeed,  suppose  the  sequence 

{u  C  D  such  that  u  -►u®  .  Well  posedness  provides  that  the 

n  1  n 

functions  s(‘)u  eC([0,'»);X)<^c'([0,'”);X)  converge  uniformly  to 
n 

S  ( •  ) u ® £ C ( [  0  ,  ” )  ;  X)  which  may  not  be  a  solution  in  the  sense  of 
Definition  2.5,  but  which  will  be  called  a  generalized  solution  if 


u°£D  (this  is  the  same  as  the  usual  notion  of  a  weak  solution  — 
see  (Fattorini,  1983:  30,31)). 

It  is  difficult,  in  general,  to  determine  whether  an  abstract 
Cauchy  problem  is  well  posed  and,  hence,  whether  there  exists  a  propa¬ 
gator  for  an  arbitrary  mathematical  model  with  the  form  of  equations 
(2.1),  (2.2).  If  the  linear  operator  A  in  equation  (2.1)  satisfies 
certain  conditions,  however,  the  propagator  can  be  shown  to  exist  and, 
iterative  schemes  are  known  for  its  construction.  Specifically,  if 
AeB(X),  or  if  A  is  in  any  of  the  G-classes  defined  previously,  then 
the  propagator  exists  and  can  be  constructed  by  an  iterative  process. 
The  details  of  this  assertion  are  now  presented. 
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whose  every  element  is  invertible.  (S^,a^)  introduced  above  is  not 
a  group  since  only  the  identity  element  of  is  invertible.  By 
enlarging  to  include  all  rational  numbers  greater  that  zero,  a 

group  can  be  constructed.  For  the  set  =  {r:r  =  p/q,p,q  =  1, 
2  ,  .  .  . }  it  is  straightforward  to  show  that  (S^  ,  )  is  a  group. 

Consider  now  the  set  5  =  {S(t);t^0}  ,  where  S  is  the 
propagator  of  a  well  posed  abstract  Cauchy  problem.  Let  a  binary 
operation  cs  be  defined  on  5xS  by 

a(S  (t,  ) ,S  (t  ) )  =  S(t  )os(t  )  (t  ,t  >0) 

where  the  symbol  "o  "  represents  the  composition  of  two  functions. 
The  following  theorem  summarizes  several  important  properties  of  the 
propagator  S  (Fattorini,  1983:  63): 

Theorem  gA. 

If  S  =  {s(t)  :t^0}  ,  where  S  ia  the  propagator  for  a  well 
posed  abstract  Cauchy  problem,  and  a  is  the  binary 
operation  defined  above,  then 

(1)  (S,a)  is  a  monoid 

(ii)  for  ,S(tj+t^)  =  S(tj)os(t^) 

(iii)  the  operator  S  :  [ 0  ,<» )^B  (X)  is  strongly 

continuous  at  every  t>0  ,  and  strongly  continuous 
from  the  right  at  t  =  0, 

(iv)  there  exist  nonnegative  constants  M,  8  such  that 

||s(t)||  < 

The  term  strongly  continuous  group  is  defined  in  a  similar 
fashion.  Consider  a  set  5'  =  fS(t):t£R}  ,  and  let  a'  be  the 


binary  operation  defined  by 


a'  (S(tj )  .SCt^))  =  S(c^  )os(t^) 

If  this  composition  satisfies 

S(tj+t^)  =  S(tj )0S(t^) 

for  all  ,t^eR  ,  and  if  S  is  strongly  continuous  at  t  =  0  ,  then 
S  is  referred  to  as  a  strongly  continuous  group  (Curtain  and 
Pritchard,  1977:  149;  Fattorini,  1983:  81). 

The  follcwing  theorem  is  the  most  important  result  in  the  study 
of  the  abstract  Cauchy  problem.  It  provides  necessary  and  sufficient 
conditions,  in  terms  of  the  operator  A  and  its  resolvent  r(z  ,  A),  for 
the  abstract  Cauchy  problem  to  be  well  posed  (Fattorini,  1983:  65). 

Theorem  2.2 

Let  the  operator  A  in  equation  (2.1)  be  closed.  The 
abstract  Cauchy  problem  (2,1),  (2.2)  is  well  posed  and  its 
propagator  S  satisfies 

||s(t)l|  £  Me®*^  (t>0) 

if  and  only  if  AeG  (M,S  )  . 

A  similar  result  for  Jie  abstract  Cauchy  problem  on  the  whole  real 
line  exists  (Fattorini,  1983:  72); 


Theorem  2.? 

Let  the  operator  A  in  the  abstract  Cauchy  problan 


u  ( t )  =  Au ( t ) 

t<oo 

(2.4) 

u(0)  = 

u°eD(A) 

(2.5) 

be  closed, 
propagator  S 

This  Cauchy  problem 
satisfies 

is  well  posed 

and  its 

||s(t)||  <  Me^l 

-“><  t  <°° 

if  and  only  if  Ac  G  '  (M  ,  3  )  . 

It  is  useful  now  to  state  a  definition  and  some  results  from 
semigroup  theory.  Theorems  2.2  and  2.3  are  usually  difficult  to  apply 
directly,  but  the  results  below  improve  the  situation  somewhat. 

Definition  2.7  (Infinitesimal  Generator)  (Curtain  and 
Pritchcird,  1977!  150,151;  Fattorini,  1983:  81) 

Let  S  =  {  S  ( t )  :  t^OlCB  (X)  be  a  strongly  continuous 
semigroup.  The  operator  A  defined  by 

,  1 im  S ( t ) u-u 

Au  =  „ ,  - 

whenever  the  limit  exists,  is  the  infinitesimal  generator  of 

5  . 


The  phrase  "  A  generates  a  strongly  continuous  semigroup  S  "  is 
frequently  used  to  mean  that  A  is  the  infinitesimal  generator  of  5, 
The  following  two  theorems  are  proven  in  (Fattorini,  1983:  81-83): 


The  linear  operator  A  generates  a  strongly  continuous 

semigroup  {s(t):t^0}  ,  with  the  property 

l|s(t)||  £ 

if  and  only  if  AsCCM.B)  . 

Theorem 

■pie  lirjear  operator  A  generates  a  strongly  continuous  group 
S  ( t )  : t<“;  ,  with  the  property 

||s(t)|l  ^  Me”' 

if  and  o  nly  if  A  c  G  '  ( M  ,  3  )  . 

Summarizing  the  results  thus  far,  it  is  apparent  that  the  problem 
of  showing  the  operator  A  in  equation  (2.1)  (equation  (2,4))  to  be  an 
element  of  C(M,3)  (  G'(M,3)  )  is  equivalent  to  showing  the 

(corresponding)  abstract  Cauchy  problem  to  be  well  posed.  In  order  to 
go  further  and  actually  solve  a  well  posed  Cauchy  problem,  one  needs 
to  construct  the  semigroup  generated  by  A  since  this  semigroup  is  the 
propagator  for  the  problem  and  provides  the  solution  u(t)  =  S(t)u° 
of  equations  (2.4)  and  (2.5),  Several  special  cases  are  now 
considered  for  the  operator  A  in  equation  (2.4), 

A B  ( X  )  .  Construction  of  the  semigroup  operator  is  most  easily 
accomplished  when  A  in  equation  (2.4)  is  an  element  of  B(X)  .  (it 
can  be  shown  that  A£B(X)  implies  that  A£:G’(M,3)  .)  The  following 
result  follows  directly  from  a  theorem  in  (Belleni-Morante,  1979: 


131): 


Theorem  2.6 

If  AeBCX),  then  A  generates  the  strongly  continuous  group 
S(t):-^<t<=“;  with  S(t)  defined  by 


S  (t) 


1  im 

n->oo 


j=o  ^  ■ 


It  can  also  be  shown  (Belleni-Morante,  1979:  130-133)  that 
satisfies 


|s  (t)l 


< 


(  -00<  t  <“  ) 


In  light  of  Theorem  2.5  and  the  foregoing,  it  is  clear  that  the 
solution  of  the  abstract  Cauchy  problem  (2.4),  (2.5)  is  given  by 

u (t)  =  S (t)u° 

for  any  u  °  £  X  . 

A:  C  ( 1  ,  0  )  ,  Ac  G  '  (1,0)  .  Consider  next  the  case  where  AeG(l,0) 
in  the  Cauchy  problem  of  equations  (2.1),  (2.2).  Define  a  sequence  of 

^  ,  OO 

operators,  by 

n 

(c>0,  n=l,2,...) 


S  (t)  = 

n 


(I--A)-' 

n 
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Ki.c)  =  ^77^7^  f  £f(x,£.t)d'p  (3.3) 

VCx.t)  =  n(x^T)  J  X^^P.)  f  (2S»£>  t )  d  S  (3.4) 

R^ 

where  v(£)  =2/7111  ,  The  current  density  vector^  ^(x,t)  ,  is  given 

by 


J(x,t)  =  qV(x,t) 

Finally,  the  pressure  tensor.  P  ,  is  defir^ed  as  follows: 


P(2<.t)  =  /  L2-P(x,t)][v(2)-V(x,t)]'^f  (x,£,t)dS  (3.5) 

The  preceding  definitions  can  aill  be  expressed  rigorously  within 
the  context  of  probability  theory.  Let  Q  =  R®  ,  and  denote  the  Borel 
field  (Maybeck,  1979s  62)  associated  with  by  Define  next  a 

set  function  for  every  t£[0,T]  by 


Pj.  ( {  (x,£)  :  (x,£)eB}  ) 


f (x ,£, t ) d 


xd 


P 


where  B-z  F  ,  and  f  is  a  distribution  function  as  defined  above.  For 
each  ts[0,T]  ,  the  triplet  (Fc,F,p^)  forms  a  probability  space 
(Maybeck,  1979:  64).  Defining  a  new  function  f*  by 

f’*(x,£,t)  =  ■|■f(x,£,t) 
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effects.  This  approach  becomes  unwieldy  for  very  large  numbers  of 
particles,  but  it  is  sometimes  taken  (Cohen  and  Killeen,  1983:  59). 

Generally,  however,  only  macroscopic  quantities  are  of  interest,  as 
opposed  to  the  specific  path  of  any  single  particle.  Consequently, 
models  of  a  plasma  usually  incorporate  probability  concepts. 

The  kinetic  theory  of  plasmas  is  frequently  developed  by  use  of  a 
distribution  function*  (Davidson,  1974:  11,  12;  Reif,  1965:  494,  495; 
Krall  and  Trivelpiece,  1973:  5,6;  Chen,  1974:  199,  200).  Suppose 

there  exists  a  collection  of  N  charged  particles  and  a  function, 

f  :  R'’ X  [  0  ,  T  ]-^  [0  ) .  f  is  called  a  distribution  function  if  the  pro¬ 
duct  f  (2<  t )  d  ^  xd  ^  p  yields  the  mean  number  of  particles  in  the 

hypercube  d  xd  p  centered  at  (x,2.^  at  time  t  . 

By  integrating  out  the  dependence  of  f  on  the  momentum  coordinates, 
the  number  density,  n  ( jc ,  t )  ,  is  obtained: 

n(x,t)  =  f f(x,£,t)d^p  (3.2) 

If  the  particles  each  have  charge  q,  then  the  charge  density, 
cr  (x  ,  t )  ,  is  given  by 


a (x , C )  =  qn  (x , t ) 


The  macroscopic  momentum  vector.  ZCii.t)  ,  and  the  macroscopic 
velocity  vector.  ,  are  defined  as  follows: 

*The  reader  is  cautioned  that  the  phrase  "distribution 
function"  in  plasma  physics  literature  is  not  synonymous 
with  a  "cumulative  distribution  function"  in  probability 
theory. 
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Notation  and  Definitions 

Most  of  the  notation  in  this  chapter  corresponds  to  that  commonly 
found  in  plasma  physics  texts.  A  summary  is  given  in  Appendix  B,  but 
for  the  reader  who  is  unfamiliar  with  this  area,  a  discussion  of  some 
of  the  pertinent  notation  and  definitions  is  now  given. 

Single  Particle.  Consider  a  particle  of  mass  ni  and  charge  q  in 
the  presence  of  an  electric  field  ^  and  magnetic  field  B  .  The  force 
on  the  particle  exerted  by  these  fields  is  given  by 

i;(t)  =  q[E.(x,t)  +  v(t)x^(x,t)] 

where  v(t)  denotes  the  velocity  vector  of  the  particle.  The  relativ¬ 
istic  version  of  Newton's  second  law  of  motion  is 

d  £(t)  =  F(t) 
d  t 

where  _p(t)  is  the  mechauiical  momentum  vector  of  the  particle.  The 
momentum  vector  is  related  to  v(t)  ^y 

£(t)  =  Ymv(t)  =  (l-6^)"‘""mv(t)  (3.1) 

where  B=lS.|=!v(t)|/c  ,  and  c  is  the  vacuum  speed  of  light. 

Plasma.  Now  consider  a  collection  of  charged  particles.  It  is 
always  possible  to  write  the  equations  of  motion  for  each  individual 
particle  including  inter- particle  forces  as  well  as  external  field 
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to  the  summary  of  notation  given  in  Appendix  B  in  lieu  of  reading  the 
following  section. 

The  most  complete  description  of  a  collisionless  plasma  consists 
of  the  self-consistent  Vlasov  and  Maxwell  equations.  Models  developed 
directly  from  these  are  known  as  microscopic  descriptions  (Davidson, 
1974:  10).  These  equations  are  presented  and  a  linear  perturbation 
model  is  developed.  This  model  is  then  shown  to  have  the  structure  of 
an  abstract  Cauchy  problem. 

By  "taking  moments"  of  the  Vlasov  equation  one  can  develop  a 
chain  of  equations  which  are  commonly  referred  to  as  macroscopic 
descriptions  (Davidson,  1974:  14).  The  continuity  and  momentum 
equations  are  the  first  and  second  set  of  equations  in  the  chain,  and 
these  are  presented  following  the  microscopic  model  discussions. 

The  microscopic  and  macroscopic  descriptions  are  stated  in  Car- 
tesian  coordinates  for  ease  of  exposition,  but  typically  a  cylindrical 
coordinate  system  is  more  practical  when  invoking  symmetry  conditions. 
Therefore,  a  coordinate  transformation  is  performed  on  the  macroscopic 
equations.  This  facilitates  development  of  a  particular  single  degree 
of  freedom  nonlinear  model.  This  model  is  linearized  about  an  appro¬ 
priately  chosen  equilibrium,  and  the  resulting  linear  model  is  also 
shown  to  have  the  abstract  Cauchy  problem  structure. 

The  final  section  is  concerned  with  a  further  simplification  of 
this  single  degree  of  freedom  model  which  isolates  certain  dominant 
dynamic  behavior. 
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of  models  with  the  abstract  Cauchy  problem  form.  Successful  develop¬ 
ment  of  such  models  would  invite  application  of  the  growing  body  of 
infinite  dimensional  modern  control  theory  to  particle  beam  dynamics 
problems. 

The  i  nter- parti  cl  e  forces  are  typically  classified  as  either  col¬ 
lective  or  collisional  forces  (Lawson,  1983:  2).  Collective  forces 
are  those  which  depend  only  upon  an  average  of  the  fields  of  many 
nei^boring  particles.  Collisional  forces,  on  the  other  hand,  depend 
upon  the  detailed  structure  of  the  charge  distribution.  The  models 
developed  in  this  chapter  deal  only  with  the  case  in  which  collective 
forces  dominate.  Collisional  forces  are  not  considered  since  particle 
accelerators  are  generally  designed  to  have  low  collisional  frequen¬ 
cies. 

The  term  "plasma"  has  been  defined  in  various  ways  in  the  litera¬ 
ture  (Lawson,  1977:  3).  In  the  present  work,  any  collection  of 
charged  particles  whose  collective  forces  are  not  ne^igible,  when 
compared  with  forces  exerted  by  external  fields,  is  termed  a  plasma. 

In  many  applications  particle  beams  are  produced  and  transported 
some  distance  in  a  vacuum.  All  models  in  this  chapter  are  developed 
under  this  assumption.  Consequently,  the  assumptions  made  thus  far 
can  be  simply  stated  as  follows:  this  chapter  is  devoted  to  the  pres¬ 
entation  and  development  of  dynamic  models  of  a  collisionless  non¬ 
neutral  plasma  in  a  vacuum. 

Overview.  Notation  and  definitions  from  electrodynamics  and 
plasma  physics  are  stated  first.  Since  the  notation  is  essentially 


standard,  the  reader  familiar  with  these  two  areas  may  wish  to  refer 


III.  Modelling  the  Dynamic  Behavior  of 
Intense  Charged  Particle  Beams 


Introduction 

Problem  Description  and  General  Assumptions.  In  this  research 
charged  particle  beam  is  considered  to  be  any  collection  of  charged 
particles  having  gross  motion  approximately  parallel  to  some  curve. 
The  curve  is  called  the  axis  and,  in  general,  the  cross-sectional 
shape  of  the  beam  varies  along  the  axis.  A  wide  range  in  complexity 
of  beam  models  exists  due  to  the  fact  that  the  particles  are  charged. 

For  a  sufficiently  low  number  density,  the  trajectories  of 
charged  particles  are  unaffected  by  the  presence  of  other  particles 
around  them.  In  this  case,  the  modelling  process  is  relatively 

strai^tf orward  since  overall  beam  behavior  can  be  inferred  from 
motions  of  individual  particles.  The  study  of  trajectories  of  parti¬ 
cles  in  low  density  particle  beams  is  referred  to  as  "charged- par  tide 
optics"  (Lawson,  1977:  3)  and  is  not  considered  here. 

Inter-partide  forces  cannot  be  ignored  at  hi^  number  densities; 
far  more  complex  and  interesting  models  are  required  in  this  case. 
Most  often  these  models  consist  of  partial  differential  equations. 
Consequently,  the  study  of  the  dynamic  behavior  of  beams  whose  inter- 
particle  forces  cannot  be  neglected  is  a  ripe  area  for  the  development 
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The  spectrum  of  A  can  be  shown  to  be  the  empty  set  in  this  case. 

The  Convection  Operator  (Belleni-Morante,  1979;  340-344). 

Consider  next  the  operator  A  :V  (A)->-X  defined  by 

Af  =  -v^f 
dx 

with  X  =  L^(-<»,o°)  and  0  (A )  =  f  f  eX  e  }  .  it  can  be  shown 

that  Ae  G '  (1,0)  and,  further,  that  the  strongly  continuous  group 
{  S  ( t  ):-“<  t<°°} generated  by  A  is  characterized  by 


S(t)u°(x)  =  u“ (x  -  vt) 


Summary 

This  chapter  has  provided  a  necessary  frame  of  reference  for  the 
next  three  chapters.  Some  notation  and  fundamental  definitions  were 
presented  first,  along  with  several  references.  Next  the  structure 
and  some  key  concepts  associated  with  the  abstract  Cauchy  problem  were 
introduced.  The  link  between  operator  semigroup  theory  and  the  Cauchy 
problem  was  then  established,  along  with  several  Important  results. 
Finally,  some  familiar  operators  were  covered  in  the  semigroup  theory 
setti  ng. 


f  a  indicated  in  the  first  chapter,  a  wide  variety  of  partial 
differential  equation  models  have  been  established  to  describe  the 
dynamic  behavior  of  a  beam  of  charged  particles.  The  next  chapter  is 


llluatrate  the  wide  applicability  of  the  theory. 

.a*~^-order  Matrix.  Let  A  =  a,  an  n*’^-order  matrix  of  real  numbers 
with  0(A)  =  R^  =  X  ,  In  this  case,  the  associated  abstract  Cauchy 


problem 


^x(t)  =  Ax(t) 


x(0)  =  x' 


with  X  =  (x^  ,  x^  ,  .  .  .  x^)  ,  is  a  finite-dimensional  model 

(ditn(X)  =  n),  and  AeB(X)  .  Consequently,  by  Theorem  2.6,  A 
generates  the  strongly  continuous  group  {  S  ( t  )'•-“<  t<«>}  where 


n 

S(t)  .  11„  y  E  e*' 


The  spectrum  of  A  consists  of  the  n  (or  fewer)  complex  numbers  X 


for  which 


det(Xl-A)  =  0 


jerator  ( Bell eni-Mor ante,  1979:  136-138).  Let  the 


operator  A  be  defined  by 


y) f (y) dy 


for  every  fex  =  C[0,l].  it  is  not  difficult  to  show  that  AeBCX) 
and  that  A  .  The  strongly  continuous  group  of  operators 
{  S  ( t )  : -°°<  t<°°}  is  defined  by 


S(t) f  =  f  + 


/  1 2  sin^ — ^ — ^^Af  +  12/^l-cos — ^ — ^A‘ 


The  operator  A  is  linear  eind  has  domain,  V  (A)  ,  in  the  Banach  space 


X  .  The  function  g  takes  on  a  value  in  X  for  each  t^O  .  The 
following  theorem  provides  sufficient  conditions  for  this  problan  to 
have  a  unique  solution  (Fattorini,  1983:  87): 

Theorem  2.12 

Let  the  operator  A  in  equation  (2.6)  be  an  element  of  the 
class  If  g  is  a  continuously  differentiable 

function  on  the  interval  [0,T]  ,  then  the  unique  solution 
of  equations  (2.6),  (2.7)  is  given  by 

u(t)  =  S(t)u'> 

where  f  S  ( t )  :  t^O } 

generated  by  A. 

Some  Familiar  Operators 

Various  operators  which  eu?e  familiar  to  engineers  and  physicists 
have  been  analyzed  in  the  literature  from  the  semigroup  operator  point 
of  view.  Results  are  now  given  for  the  following  operators:  (1)  an 
n*''^-order  matrix  of  real  numbers,  (2)  a  specific  integral  operator, 
and  (3)  the  (scalar)  convection  operator.  The  first  example  involves 
a  bounded  operator  defined  on  a  finite  dimensional  space,  the  second  a 
bounded  operator  defined  on  an  Infinite  dimensional  space,  and  the 
final  example  deals  with  an  unbounded  operator  defined  on  an  infinite 
dimensional  space.  These  specific  operators  are  chosen  solely  to 


+  j S  ( t- s  )  g  ( s ) d s  (0£t£T) 

0 

is  the  strongly  continuous  semigroup 


following  result  is  very  practical  (Belleni-Morante,  1979!  179-182): 


Theorem  2.11 

If  A  =  A  +A  A  eB(X)  ,  and  A  eC(M,3)  then 
AeC(M,3+M  A^  ) . 

D 


As  a  resialt  of  this  theorem,  one  need  "worry"  only  about  the 
"unbounded  portion"  of  an  operator,  usually  the  derivative  terms.  The 
strongly  continuous  semigroup  for  the  operator  A  satisfying  this 
theoran  is  constructed  by  an  iterative  process.  Let  {s(t):t^0}  be 
the  semigroup  generated  by  A^  ,  and  define  the  sequence 


(t)f  =  S(t)f  (t>0) 


Z^^^(t)f  =  S(t)f  +  /S  (t-s)Aj^Z^(s)  f  (s)ds  (t>0,  n=l,2,...) 


The  strongly  continuous  semigroup  {z(t):t^0}  generated  by  A  ,  then, 
is  defined  by 


Z  (t ) f  =  lira  Z . (t) f 

j-voo  J 


(VfeX,  t>0) 


The  final  result  in  this  section  is  a  theoran  concerning  the 


inhomogeneous  problem: 


=  Au(t)  +  g(t) 


(2.6) 


u(0)  = 


(2.7) 
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(2.4),  (2.5))  is  given  by 


u(t)  =  Z(t)u° 

for  any  u°£p(A)  and  for  all  t^O  (-oo<t<oo 


Further  Practical  Results 

Three  results  which  are  often  of  use  in  the  application  of  the 
theory  are  now  introduced.  The  first  is  generally  useful  if  the 
underlying  space  X  is  a  Hilbert  space  and  the  norm  corresponds  to  the 
energy  of  the  system. 

Theorem  ZASl 

Let  A:P(A)^X  ,  where  P(A)CX  and  P(A)  is  dense  in  the 
Hilbert  space  X  .  Then  A£G(1,3)  if  and  only  if 

(i)  (zI-A)P(A)=X  Vz  such  that  Re(z)>B 

(ii)  Re(Af  ,  f)<8  Ilf  11^  VfeP(A) 

A  densely  defined  linear  operator  A  satisfying  condition  (ii)  with 
6  =  0  is  called  dissipative!  also,  if  “A  is  dissipative  then  A  is 
cauLled  accretive.  For  further  discussion  in  this  area  and  a  proof  of 
Theorem  2.10,  refer  to  ( Bell eni-Mor ante,  1979:  142-145). 

Often  a  complicated  operator  A  can  be  broken  into  two  opeators: 
A  =  A  +A  .  If  A  is  chosen  such  that  it  is  a  bounded  linear 

DU  D 

operator,  defined  on  all  of  the  underlying  Banach  space  X,  then  the 
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A e  G  ( M  ,  3  )  .  Ae  C  '  (M  ,  3  •).  The  proofs  of  Theorems  2.7  and  2.8  are 


easily  modified  for  AeG(M,0)  orAeG'CM.O)  .  Furthermore,  if 
AeG(M,3)  (AeC  (M,3)  )  then  the  operator  A^  ,  defined  by 

A^  =  A-3I 

is  in  the  class  G(M,0)  (  G'  (M,0)  ).  Consequently,  the  following 
theorem  can  be  proven  with  little  additional  work  ( Bell eni-Mor ante, 
1979:  159): 

Theorem  2.Q 

If  A£G(M,S)  ,  then  A  generates  the  strongly  continuous 
semigroup  {z(t):t^0}  with  Z(t)  defined  by 

6  t 

Z(t)u  =  S(t)u  (VueX,  t^O) 

where  {S(t):t^0}  is  the  semigroup  generated  by 

A^  =  A-3I. 

The  ar^logous  result  for  AeG'(M,B)  follows  immediately.  The  norm  of 
Z(t)  satisfies 

llz(t)lj  ^  (-»<[<'») 

for  AeG(M,S)  ,  and,  if  AeG'(M,3)  generates  the  strongly  continuous 
group  f  Z  ( t )  :  -“<  t  <“’ }  , 

11  Z  ( t )  II  £  (_c»<t<oo) 

In  either  case,  the  solution  of  equations  (2.1),  (2.2)  (equations 
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Theorem  2.7 

If  A£C(1,0)  ,  then  A  generates  the  strongly  continuous 

semigroup  with  s(t)  defined  by 


S(t)u  =  (VueX,  -oo<t<'») 


Additionally,  this  semigroup  satisfies 


|S(t)|l  <  1 


and,  hence,  the  solution  to  the  Cauchy  problem  of  equations  (2,1), 
(2.2;  is  again 


u(t)  =  S(t)u“  (t>0) 

for  any  u’  sP  (A  )  . 

Letting  S^(c)  be  defined  as  above,  but  for  -oo<t<oOj  one  also 
has  the  following  ( Bell eni-Mor ante,  1979!  160): 


Theorem  2.8 

If  AeG'CljO)  ^  then  A  generates  the  strongly  continuous 
group  { S  ( t  ):-«><  t<'”}  with  S(t)  defined  by 

S(t)u  =  lim  S^(t)u  (VueX,  -oo<t<co) 

The  group  thus  defined  satisfies 

l|S(t)||  <  I 

and  the  solution  of  equations  (2,4),  (2,5)  is 

u(t)  =  S(t)u'> 


for  any  u  °  e  P  (A )  , 
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it  is  straightforward  to  show  that  f*  is  a  joint  probability  density 
function  (pdf)  for  each  te[0,T]  with  random  variables  k  and  £  , 
Integration  of  f*  over  all  £eR^  yields  a  marginal  probability  den¬ 
sity  function; 


=  /  f*(x,£,t)d^p 

This  marginal  pdf  is  related  to  the  number  density,  n(x,t)  ,  by 

n(x,t)  =  Nf*(x,t) 

X 

A  conditional  probability  density  function  is  now  needed  to  express 
the  macroscopic  momentum  vector,  the  macroscopic  velocity  vector,  and 
the  pressure  tensor  in  terms  of  the  probability  space.  Specif  iced,  ly, 
let  the  function  be  defined  by 

f  * (x ,£,  t )  f (x ,£,  t ) 

f*i  (£;x,t)  =  ITT - ^ ^ - 

£|X  ^  -  f*(x,t)  n(x,t) 

Since  ^£|  X  ^  every  ( x  ,  t )  e  R^  x  [  o  ,  t]  ,  conditional 

expected  values  of  any  function  6(p)  can  be  taken; 

E[''r(p)[x;  =  j  ':'(p)fp|^(p:x.t)  d'p 

R3 

Conditional  expected  values  of  the  functions  £  ,  ^(P.)  >  and 

[£  -  £(2C,t)][v(£)  -  VCx.t)]"^ 

yield  P(£,t)  ,  respectively. 
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Expression  of  the  basic  definitions  of  plasma  physics  in  a  proba¬ 
bility  theory  setting  provides  rigor  and  clarity  for  applied  mathema¬ 
ticians.  On  the  other  hand,  the  notation  used  by  plasma  physicists  is 
both  intuitive  and  well-established.  Consequently,  now  that  the  con¬ 
nection  between  these  two  areas  has  been  established,  plasma  physics 
notation  and  definitions  are  used  in  the  ranainder  of  this  disserta¬ 
tion. 

The  rationalized  MKS  system  of  units  is  used  in  this  chapter 
since  this  seems  to  be  the  choice  of  many  authors  of  chau'ged  particle 
beam  texts  ((Lawson,  1977),  for  example),  Hcwever,  it  should  be  noted 
that  most  plasma  physics  authors  prefer  the  cgs  Gaussian  system  (for 
example,  see  (Davidson,  1974;  Krall  and  Trivel piece,  1973)).  Both 
systems  have  their  advantges  and  disadvantages,  and  the  transition 
from  one  system  to  the  other  is  not  difficult.  In  the  rationalized 
MKS  system,  the  symbols  £  ,  p  are  used  to  represent  the  absolute 
dielectric  constant  and  magnetic  permeability  which  are  related  by 


Finally,  since  vector  cross  products  are  somewhat  tedious  to 
write  out  in  detaiil,  the  permutation  symbol  (Marion  and  Heald,  1980: 
456),  z  .  as  defined  in  Appendix  B,  and  summation  notation  are  fre- 
quently  used.  By  way  of  example,  consider  the  cross  product 

W  =  V 

The  components  of  w  can  be  expresed  compactly  as 
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w  . 
1 


G  .  .  ,  U  .  V, 
ijk  j  k 


(1=1.2, 3) 


For  instance,  if  i=l,  the  above  expression  yields 

3  3 

"'l  ^  Z  Z^ljk  =  -  ^^2 

j=l k=l 

since  £,,=1,  e,,=-l  ,  and  e ,  =0  for  all  other  possible  trip- 

i23  132  Ijk 

les  ( 1  ,  j  ,  k). 


Microscopic  Descriptions 

Vlasov  Equation.  The  distribution  function  of  a  nonneutral  col- 
lisionless  plasma  of  a  sin^e  species  obeys  the  Vlasov  equation 
(Davidson,  1974:  12).  If  ^(x.t)  represent  the  total  elec¬ 

tric  and  magnetic  field  at  time  t,  the  Vlasov  equation  in  Cartesian 
coordinates  can  be  written  as 

+  q[E.(x,t)  +  E  .  .  B  (x  ,  t  )  (x  t  )  =  0  (3.6) 

The  fields  £  ,  _B  arise  from  external  charges  as  well  eis  from  collec¬ 
tive  effect  from  the  particles  in  the  plasma  itself.  Denoting  the 
external  fields  by  and  the  self  fields  by  »  the  total 


fields  can  be  expressed  as 


Maxwell's  Equations.  Maxwell's  equations  must  be  satisfied  as 
well  as  the  Vlasov  equation.  The  external  fields  are  produced  by 
external  charges  or  current  densities,  but  since  these  will  ultimately 
be  regarded  as  controls  which  can  be  applied  in  a  prescribed  manner, 
their  corresponding  Maxwell  equations  are  unimportant  at  present.  On 
the  other  hand,  the  self  fields  depend  intimately  upon  the  distribu¬ 
tion  function  through  Maxwell's  equations: 


3c  ^ 


3  X  . 


(3.7) 


9  B  .  (  X  ,  t )  = 


j 


(3.8) 


3 _ E"(x,t) 

9x. 


=  a (x, t)  /£, 


(3.9) 


9 _ B  (x,t) 

9x  ^ 

J 


(3.10) 


for  i-=l,  2,3.  Recalling  that  o  and  ^  depend  upon  the  distribu¬ 
tion  function  f  ,  it  is  seen  that  (3.6)  through  (3.10)  represent  a 
system  of  nine  coupled  nonlinear  integro-diff erential  equations.  Rep- 
resenting  the  ordered  set  (f(x.£.t),  ^  (x,t),  ^  (^c.t))  by 
u(t)  ,  equations  (3.6)  through  (3.8)  can  be  written  as 


d  u(t)  =  F(E®,B®) (u(t) ) 


(3.11) 
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where  F.  depends  upon  the  external  fields  and  represents  the  nonlir»ar 
operations  indicated  in  those  equations.  Furthermore,  equations 

(3.9) »  (3.10)  serve  as  restrictions  on  the  domain  of  F. ,  as  would 
boundary  conditions  which  are  typically  present  in  any  given  physical 
situation,  A  solution  of  the  differential  equation  (3.11),  and  an 
associated  initial  value,  u_(0)  =  u.**  ,  is  generally  difficult  to 
obtain. 

Linearization  of  the  System.  If  the  nonlinear  operator  is 
approximated  by  a  linear  operator,  the  resulting  system  can  be  shown 
to  be  an  abstract  Cauchy  problan.  This  is  now  demonstrated  for  the 
special  case  of  both  the  the  electric  and  magnetic  external  fields 
being  identically  zero. 

In  infinite-dimensional  systans,  nonlinear  operators  can  be 
approximated  in  a  manner  analogous  to  the  first-order  Taylor  series 
technique  in  finite-dimensional  systems.  Consider  the  equation 

x(t)  =  ^(x(t)) 

where  g  is  a  vector-valued  nonlinear  function;  if  jc®  is 

known  to  be  a  solution  of  £.(2L'’)=0f  and  if 


then  x°  is  called  an  equilibrium  solution.  Suppose  ^  can  be 
represented  by  a  Taylcr  series  at  jc”  ; 


^(x)  =  _£(x'’  )  +  )  (•i~2i'’  )  + 
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the  original  system  can  be  approxi- 


Now,  letting  ^(t)=x(t)-x°  , 
mated  by  the  linearized  perturbation  equation 

^(t)  J^oAx(t)  (3.12) 

where  J  o  =  ^  S.Cii'’).  Th®  operator  J  ois  the  Frechet  derivative  of 

—  Sjc,  — 

the  nonlinear  £  at  x°  provided  each  entry  in  J  o  is  continuous  (see 

X 

examples  1  and  4  of  (Luenberger,  1969:  171-174)).  In  li^t  of  the 
comments  in  Chapter  II  following  the  discussion  on  Frechet  deriva¬ 
tives,  the  Gateaux  derivative  of  is  also  J^o  • 

In  many  situations  it  is  not  possible  to  show  that  an  operator  is 
Frechet  differentiable,  although  the  Gateaux  derivative  can  usually  be 
determined  in  a  straightf orwau*d  manner.  Linearizations  based  on  the 
Gateaux  derivative  cannot  be  justified  rigorously  a  priori,  but  such 
models  are  often  used.  Solutions  obtained  for  these  models  should  be 
verified,  if  possible,  by  alternate  methods. 

Let  the  operator  ,  X  a  Banach  space,  represent  the 

Gateaux  derivative  of  F:D(F)-»-X  at  u”,  and  consider  the  equation 

^u(t)  =  F(u(t))  (3.13) 

d  t 

Suppose  u°  is  an  equilibrium  solution  (defined  in  the  same  manner  as 
in  the  finite-dimensional  case:  F(u°)  =  0)  and  let  6 u  ( t )  =u  ( t ) -u°  . 

Provided  the  Gateaux  derivative  is  a  lineau?  operator,  the  approxima¬ 
tion 
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l3  made 


This  is  a  direct  result  of  the  definition,  since  if 
5u(t)=hv(t)  ,  for  any  veX,  heR,  the  following  limits  exist: 

lim  F  (u°+hv  (t ) )-F  (u° )  =  lim  F ( u “ +6 u ( t ) ) - F (u °  )  =  1 im  15F  nSuCt) 
h-0  h  h -  h-vO  u 

By  definition  of  the  limit,  then,  this  implies 


Consequently,  for  h  sufficiently  small,  the  differential  equation 

d_5u(t)  =  5F  o'SuCt) 
dt 

becomes  the  infinite-dimensional  analogue  of  the  finite-dimensional 
result  of  equation  (3.12), 

With  the  approximation  (3.14),  the  nonlinear  Vlasov  equation 
(3.6)  can  be  linearized  in  a  straightforward  manner.  Let 

,  b"  be  equilibrium  solutions  of  equations  (3.6) 
through  (3.10),  and  define  5f,  §_E,  6_B  in  the  obvious  way.  Then  the 
lineau'  approximation  to  the  Vlasov  equation  is 


3 _ 6f(2<,2_,t)  +  V  .  3 _ <Sf(x,£^,t) 

3t  ^3x. 


+  q(E°(x)  +  e  .  . ,  V  .  b“  (x  ,  t  )  )3 _ 6f(x,£.,t) 

1  ijk  3  k  -  g 

1 

+  q  (5  E  .  (2i,£,  t )  4-  e  .  . ,  V  .  6  B  (x  ,  t  )  3 _ f“(x,£)  =  0 

1  IJ  K  3  K  g 

i 


(3.15) 


Since  Maxwell's  equations  are  linear,  simply  replacing  the  functions 
(  f  ,  E^  ,  )  with  (5f  ,  in  equations  (3.7)  through  (3.10) 

yields  the  linearized  versions  of  ty;ese  equations: 


3 _ 5  E  .  (jc,  t ) 

3t  ^ 


=  f  v^5f(x,£,t)d^p+c^^  3 _ 6B®(x,t)  (3.16) 

J  ,  3  X  . 


^6B®(x,c)  =  -£..,3 _ E^(x,t) 

3t  ^ 


3 _ 5  E  .  (x  ,  t )  =  q 

3  X  .  e  „ 

J  “ 


3 _ 5B“(x,t)  =  0 

3  X  . 

J 


J'sf  (x  ,£, 


t )  d  p 


(3.17) 


(3.18) 


(3.19) 


The  system  of  equations  (3.15),  (3.16)  and  (3.17)  represents  an 
abstract  Cauchy  problem 


d  w(c)  =  Aw(t)  (t>0) 

d  t 

with  initial  condition  w(0)  =  w ” .  The  underlying  Banach  space 

7 

X  =  u  X .  is  yet  to  be  specified.  The  domain  of  A  should  include 

i=  1 

the  restrictions  of  equations  (3.18),  (3.19),  as  well  as  any  addi¬ 
tional  boundary  conditions. 

Further  analysis  of  the  microscopic  equations  requires  realistic 


»*v 


equilibrium  solutions  which  are  smooth  enough  for  their  derivatives  in 
equation  (3.15)  to  exist.  Such  solutions  are  not  known  at  this  time, 
although  equilibrium  solutions  involving  the  Dirac  delta  function  have 
been  discovered  (e.g.,  see  (Hammer  and  Rostoker,  1970:  1831-1834)). 
Various  attempts  were  made  to  continue  analysis  of  microscopic  models 
using  such  equilibrium  solutions,  but  the  resulting  linearized  equa¬ 
tions  were  intractable. 


Macroscopic  Descriptions 

Equations  can  be  developed  from  the  Vlasov  equation  which 
describe  the  evolution  of  certain  "averaged"  quantities.  Such  equa¬ 
tions  are  termed  macroscopic  descriptions,  and  the  first  two  sets  of 
these  are  presented  below.  These  descriptions  are  appealing  since  the 
unknown  functions  associated  with  them  are  more  intuitive  than  the 
distribution  function  in  that  the  physical  quantities  involved  are 
more  directly  observable.  However,  certain  phenomena,  such  as  Landau 
damping,  cannot  be  predicted  by  such  descriptions  (Davidson,  1974: 
11),  and,  consequently  some  information  is  forever  lost  once  micro¬ 
scopic  descriptions  are  abandoned. 

The  macroscopic  equations  are  derived  by  multiplying  the  Vlasov 
equation  by  an  appropriate  function  and  integrating  over  all  momentum 
space.  Details  are  not  presented  here  since  they  can  be  found  in  var¬ 
ious  plasma  physics  texts  (see,  for  example,  (Chen,  1974:  211-213)) 
The  first  two  sets  of  equations  are  commonly  referred  to  as  the  con¬ 
tinuity  equation  (3.20)  and  the  momentum  equations  (3.21)  (Krall  and 
Trivelpiece,  1973:  88;  Davidson,  1974:  14): 


9  n (x  ,  t ) 
3  t 


(3.20) 


=  -  3  [n(x.t)V(x.t)l 

3x  . 

1 

9  P  .(ji.t)  +  V.(x,t )  3  P  .  (x  ,  t )  +  1  3  [  P  (x  ,  t )  ]  .  . 

3t  ^  3x.^  n(x,t)3x. 

1  —  1 

=  q[Ej(x,t)  +  C^j^^Vj^(x,t)B^  (x,t)  ]  (3.21) 

where  j  =  1,2,3  ,  and  [P(x,t)]^j  is  the  ( i  ,  j  )- component  of  the 

pressure  tensor  (see  equation  (3.5)). 

Equations  (3.20)  and  (3.21)  cannot  be  solved  without  knowledge  of 
the  pressure  tensor,  P  .  The  components  of  P  would  appear  as  time 
derivatives  in  the  next  higher  moment  equation,  the  energy  equation, 
but  a  quantity  would  be  needed  in  this  equation  from  the  next  higher 
moment  equation,  and  so  forth.  This  chain  of  mcanent  equations  is  fre¬ 
quently  broken  here  by  making  some  approximation  to  P  .  If  the  spread 
in  the  momentum  is  small  at  every  point,  then  components  of  P  are 
small  and,  in  the  limiting  case  of  the  momentum  being  a  deterministic 
quantity  everywhere,  P  =  0  (Davidson,  1974:  16).  The  spread  in 
velocity  is  also  zero,  in  this  case,  and  thus  the  temperature  vanishes 
everywhere.  This  idealized  case  is  termed  the  cold  plasma 
approximation. 

As  mentioned  previously,  the  macroscopic  equations  of  (3.20)  and 
(3.21)  are  somewhat  more  intuitive  them  the  Vlasov  equation  which  they 
replace.  Integral  operators  are  required  in  microscopic  descriptions 
(see  equations  (3.7),  (3.9)),  but  are  unnecessary  in  the  macroscopic 
descriptions.  Furthermore,  for  the  cold  plasma  approximation. 
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although  macroscopic  descriptious  replace  a  single  unknown  function 
(the  distribution  function,  f  )  with  four  unknown  functions  (the  num¬ 


ber  density  and  the  three  components  of  the  macroscopic  momentum  vec¬ 
tor)  the  reduction  in  the  number  of  independent  variables  is  three 
(from  (2<,_p,t)  to(2<,t)).  For  many  applications,  only  these  macro¬ 
scopic  functions  are  of  interest.  In  light  of  these  observations,  it 
is  concluded  that  macroscopic  models  are  preferred  in  the  design  of 
particle  beam  control  components,  so  long  as  they  accurately  describe 
the  number  density  and  the  macroscopic  momentum. 


A  Singlg  Desrgg  Freedom  Linear  Model 

Introduction.  Various  additional  assumptions  are  introduced  in 
this  section  in  order  to  derive  a  suitable  model  for  subsequent  illus¬ 
tration  of  semigroup  theory  techniques.  A  broad  variation  in  operat¬ 
ing  conditions  exists  for  charged  particle  beams.  Each  assumption 
below  has  been  invoked  in  plasma  physics  research  in  the  investigation 
of  beam  behavior  under  a  specific  operating  condition  (e.g.,  see 
(Davidson,  1974)).  Some  assumptions,  for  example  the  nonrelativistic 
velocity  assumption,  serve  only  to  call  out  a  specific  regime  of  oper¬ 
ation.  Other  assumptions,  such  as  the  assumption  of  the  adequacy  of 
macroscopic  descriptions,  are  made  to  simplify  the  model,  with  the 
justification  being  that  they  have  previously  been  invoked  by  plasma 
physics  researchers  and  have  been  found  to  be  useful  and  adequate  in 
describing  beam  dynamic  phenomena.  In  either  case,  the  philosophy 
taken  now  is  that  simple,  though  less  accurate  models  whose  dominant 
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u 

S 

behavior  can  be  expressed  analytically,  are  superior  in  preliminary  ’ 

control  designs  to  more  precise  models  which  require  computer  gener- 

* 

ated  numerical  solutions.  ' 

■a 

The  cold  plasma  assumption  is  in  keeping  with  this  philosophy.  As 

previously  stated,  this  is  equivalent  to  assuming  that  the  momentum  is 

deterministic  at  every  point.  In  practice,  if  the  momentum  spread  at  ‘ 

every  point  is  sufficiently  small,  then  the  cold  plasma  assumption  is 

reasonable.  Otherwise,  approximations  of  the  pressure  tensor  might  be 

required  (see  (Krall  and  Trivelpiece,  1973:  89)).  • 
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The  assumptions  made  in  no  way  limit  the  applicability  of 
semigroup  theory.  For  example,  the  same  techniques  used 
below  to  analyze  the  single  degree-of-freedom  model  could  be 
applied  to  a  three  degree-of-freedom  model.  See  Appendix  D 
for  the  structure  of  such  a  model. 
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Reduction  in  the  number  of  degrees  of  freedom  is  also  useful  in 
simplifying  the  model.  However,  models  which  allow  only  a  single  spa¬ 
tial  degree  of  freedom  in  a  Cartesian  coordinate  system  are  generally 
unrealistic,  so  a  cylindrical  coordinate  system  is  introduced. 

Since  most  applications  of  a  charged  particle  beam  require  only 
that  the  beam  operate  near  some  design  equilibrium  solution,  a  linear 
model  about  an  equilibrium  point  should  be  adequate  for  control  pur¬ 
poses.  Indeed,  if  the  control  function  desired  is  that  of  regulation, 
deviations  from  the  equilibrium  will  be  constrained  to  be  small  by  the 
action  of  the  regulator.  This  notion  is  fundamental  to  control  theory 
design  and  has  been  applied  with  success  routinely  over  the  years. 

Assumptions .  The  following  assumptions  are  in  effect  in  the 
development  below: 

A1 .  the  beam  is  in  a  vacuum 

A2.  all  velocities  are  nonrelativistic 

A3,  the  momentum  spread  at  each  point  is  small  (cold  plas¬ 
ma) 

A4.  macroscopic  descriptions  are  adequate 

A5.  the  beam  is  uniform  in  the  azimuthal  direction 

A6 .  the  beam  is  uniform  in  the  axial  direction 

A7 .  deviations  from  the  equilibrium  solution  are  small 

Additional  assumptions  are  needed  later  (page  III-24)  for  the  develop¬ 
ment  of  a  specific  equilibrium  solution,  and  are  stated  at  the  outset 
of  that  development. 

Assumptions  (Al)  through  (A4)  result  in  the  following  system  of 
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i 


can  be  computed  as  follows: 


at  z“£iJxX  ,  5g^o  , 


1  im 
h“^  0 


G_(^°  +h£) 
h 


G(z“  ) 


q 

m 


+  ~  +  z®^  -  Cz'’  +  z°C  ) 

T  Zg-,  5  -r  Z.^g  VZjj|,g  -r 


+  Z°  ;  +  z\ 

8  5  5  8 


(z“  C  +  z  t  ) 

9  4  4  9 


0 

0 

0 

0 

0 


This  expression  represents  the  Gateaux  derivative  of 
equilibrium  £  and  associated  external  equilibrium 
j?  ’  .  For  the  rigid  rotor  equilibrium,  the  vector  £ 


(3.66) 


G  for  a  general 

_  e  ,  0 

fields,  ^  , 

is  given  by 
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0 


E®(r,t)  +  B®(r,t)lJ3  - 

E^^Cr.t)  +  B®(r,t)y^  -  B^Cr.t)^^ 


G(E®,  B®)(U) 


m. 


E^(r,t)  +  B^'=(r,t)u 


B^(r,t)y3 


(3.65) 


0  0 

Since  the  external  fields  ^  are  being  viewed  as  the  controls, 

define  the  function  space  consisting  of  all  physically  attainable  con- 

6 

trols  as  y  -  u  Y.  where  E^eY,  ,  E^eY,  ,  ...B®eY,  .  Conse- 

r  ^  ’  z  ° 

0  0 

quently,  the  mapping  G(^  ,  ^  )  represents  a  unique  nonlinear  opera- 
tor  for  every  (E  ,  ^  )Ey  or,  alternatively,  one  can  consider  the 
mapping  £  as  one  which  takes  elements  from  yxX  into  X  .  This  lat¬ 
ter  interpretation  allows  an  approximation  of  ^  by  the  Gateaux  deriv¬ 
ative  in  a  manner  analogous  to  the  approximation  above  for  F  . 

Consider  an  arbitrary  element  £  in  y^-X  ,  and  suppose  the  X^, 
i  =  1  ,  2  ...  9  ,  are  all  Banach  spaces.  The  Gateaux  derivative  of  G_ 
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1  im 

u“ h-0 


F(u°+hn)  -  F(u°) 


>SF..o(f)  = 


(-  ^(ru")  + 

rdr  2  2dri 


4 


+  3.^  _  3.u°],  ,  „  , 

m  ‘  5  m  “  8  m  3  '  9 


+  ^u°n 


(—  +  ~'3q  )n2 

rdr  3  m  8  ^  '  2 


/  0  0  \ 

(— U  ^  -  -; - U  ,  )  n 

m  8  dr  8  2  12 


<7^°  “277>^3  ^  he  -  ^^9 


0  d  ,  Q  1  Q  0 

^277^^4  +  +  7^^^8 


-Uoqc'(u°n,  +  u!;n2) 


■u^qc^  (u®ni  +  ujrij)  -  c'fpn, 


-u^qc^  (u°n,  +  +  f 


c  ^  d 

-  T^rri 


dr' 


dr  2 


1  d 


r  dr '■'^6 


(3 


This  expression  is  the  Gateaux  derivative  for  F  for  a  general  equi¬ 
librium  solution  u^  ,  For  the  rigid  rotor  equilibrium  (equation 
(3.60)). 


.;-l  V'-  '."V  n 


.63) 
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merits  of  ^(t)  inequations  (3.41 )  ,(3.^2) ,  and  (3.43)  which  involve 
the  external  fields,  and  it  is  discussed  in  more  detciil  presently.  A 
linear  model  which  approximates  this  system  is  now  developed. 

Let  continue  to  denote  the  rigid  rotor  equilibrium  solu¬ 

tion  (equation  (3.60)),  and  the  z  -component  of  the  external,  uniform 
magnetic  field,  respectively.  Define  the  perturbed  variables  6  u  ( t ) 

in  the  usual  way;  =  ij(t)  -  .  Similarly,  let  the  per- 

^0  0 

turbed  external  fields,  (r,t)  ,  (r,t)  be  given  by 

^®(r,t)  =  E®(r,t)  -  E®’“(r)  =  E  ®  ( r  ,  t  ) 

6_B®(r,t)  =  B®(r,t)  -  B®’®  (r)  = 


If  all  the  spaces  i  =  1,2, ...9,  are  Banach  spaces,  then  the 

nonlinear  operator  £  can  be  approximated  by 

F(u'’  +  iu(t))  -  Ku®)  +  ^^o(5u(t))  (3.62) 

where  the  (linear)  Gateaux  derivative  of  £  at  (see  the 

discussion  associated  with  equation  (3.14)  in  the  section  "Microscopic 
Descriptions").  Calculation  of  is  straightforward: 


Bj(r  ,t) 

Bg®  ( r  ,  t ) 
B®(r,t)  -  b“ 
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1973!  117).  Finally,  for  ,  two  equilibria  aire  possible 

G  P 

(termed  the  slow  mode/  fast  mode  (Davidson,  197^!  7)). 

The  complete  solution  for  re(0,R]  is  now  summau'ized.  Letting 

u“  =  [n°(r)  ,  V“(r)  ,  Vg“(r)  ,  V®(r)  ,  E°(r)  ,  E^°(r)  , 

T 

e“  (r)  j  Bg  (r)  ^  B°  (r)]  ,  the  rigid  rotor  equilibrium  solution 


IS 


= 


n 

0 

-ojr 

V° 


iEL 


2e 


(3.60) 


2 
0 

Linearization.  Equations  (3.40)  through  (3.46)  and  (3.48) 
through  (3.50)  are  seen  to  form  a  systan  of  nine  coupled  nonlinear 
differential  equations  and  can  be  expressed  as 


^u(t)  =  F(u(t))  +  G(E®,  B®)(u(t))  (3.61) 

9 

where  F:P(F)CX->-X  ,'^=  HX  ^  and  ^(t)  =  [n(r,t),  V  (r,t), 

i=l 

Vg  (r  ,  t )  ,  (r  ,  t )  ,  (r  ,  t  )  ,  Eg  (r  ,  t )  ,  E^  (r  ,  t )  ,  Bg  ( r  ,  t )  , 
^  T  0  G 

B^(r,t;)]  ,  The  operator  £(^  ,  ^  )  represents  operations  on  ele- 
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The  term  qB°/m  ,  commonly  called  the  cyclotron  freauencv  (Lawson, 


1977!  17) »  appears  often  in  plasma  physics  and  is  denoted  byco^.  The 

term  within  the  brackets  is  simplified  by  recalling  that 

.2  .  _J_ 


0  0 


and  thus,  this  term  reduces  to 


[1  ~  (0“)^ 


By  the  nonrelativistic  assumption,(3°)  is  ne^ected.  The  phrase 
Plasma  freauencv  is  given  to  the  expression  (Lawson,  1977:  119) 


1,1 

and  it  is  given  the  symbol  <^p  .  The  following  expression  for  Vg°  (r) 
emerges  in  light  of  the  foregoing: 


( r  )  = 


_  Ir 


[o)  ±  [co^  -  203^  ] 


(3.59) 


The  use  of  the  phrase  "rigid  rotor  equilibrium"  is  justified  by  this 
expression. 

The  expression  within  the  radical  in  equation  (3.59)  provides  a 
minimum  value  for  b”.  The  external  magnetic  field  is  "sufficiently 
large"  (see  assumption  El  above)  if 

OJ'  >  >  B”  > 

C  =  P  =  Eq  -* 

0) 

?  c 

If  (jj  =  2a)_  ,  then  the  rotation  rate  is -r-,  and  the  phrase  "Bril- 

louin  flew"  is  used  to  describe  this  situation  (Krall  and  Trivelpiece, 
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^q(v;)^r 


0<  r<R 


(r) 

dr  r 


5J[L_r 

En 


0<r<R 


The  solutions  of  these  equations  are 


/  Poqn“v' 


Uoqti''v“R 


z  R 


0<r<R 


E"’“(r)  = 
r 


in^R  R 


2£o  r 


0  r  R 


Applying  these  solutions  to  equation  (3,52)  for  re(0,R] 


obtai ns 


7^ - =  0 

2m 


Solving  this  quadratic  equation  for  (r)  yields 


qO  iir^nolz  .  L2_o  U,q^(v“)^n 
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must  satisfy) . 

The  filial  assumption  simplifies  the  solution  considerably  and  is 
realistic  so  long  as  the  beam  is  nonrelativistic  (Krall  and  Trivel- 
piece,  1973s  117). 

Assumptions  (El)  through  (E5)  are  now  applied  to  equations  (3.40) 
through  (3.51).  First,  note  that  equations  (3.40),  (3.42),  (3.43) f 

(3.44),  (3.47) f  and  (3.51)  are  all  trivially  satisfied.  Letting  the 

superscript  "on  denote  an  equilibrium  solution  function,  the  remaining 
equations  become 


+  (r)+vi'  (r)B^’“  (r)-V°  (r)Bg®’° 

m  r  0  2  2  0 

(r)]  =  0 

(3.52) 

=  0 

(3.53) 

E^‘“(r) 

=  0 

(3.54) 

U  qc^V° (r)  +  c^d  B®’° (r) 

0  ^  z  —  z 

dr 

=  0 

(3.55) 

-M^qc^V°(r)  +  c  ^  d  (  r  Bg^  *  °  ( r )  ) 
r  dr 

=  0 

(3.56) 

1  l_(rE^’'’  (r)) 

=  q  n° (r) 

(3.57) 

r  dr 


Denote  the  constant  external  magnetic  field,  number  density,  and 
z  component  of  velocity  by  ,  andV®,  respectively.  The  last 

two  equations  become 
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Consider  a  beam  In  the 


shape  of  a  long  cylinder  of  circular  cross-section  with  radius  R,  the 
axis  of  which  coincides  with  the  z  axis  of  a  cylindrical  coordinate 
system.  Certain  assumptions  are  required,  in  addition  to  those  previ¬ 
ously  stated,  for  the  rigid  rotor  equilibrium  solution: 

El.  a  (sufficiently  large)  uniform  magnetic  field  in  the  z 
direction  is  the  only  external  field 

E2.  the  niMber  density  is  constant  for  and  vaxi- 

ishes  on 

E3.  the  velocity  in  the  z  direction  is  constant  for 

0^r<  R 

E4.  the  radial  velocity  is  Identically  zero 

E5.  the  z  component  of  the  self  magnetic  field  is  negligi¬ 
ble  compared  to  the  external  field 

Some  discussion  of  these  additional  assumptions  is  now  in  order. 

The  external  field  in  assumption  (El)  is  necessary  to  offset  the 
repulsive  forces  which  would  cause  the  beam  to  expand  indefinitely. 
The  particles  undergo  a  helical  motion  in  the  presence  of  this  magnet¬ 
ic  field,  resulting  in  a  balance  between  the  repulsive  forces  (elec¬ 
trostatic  and  centrifugal),  and  the  constrictive  force  (magnetic) 
whenever  the  external  magnetic  field  is  sufficiently  large.  Other 
means  of  confining  a  beam  to  a  finite  radius  are  possible  (such  as  by 
neutralization  by  background  ions  (Lawson,  19T7J  258),  fca*  example). 

Assumptions  (E2),  (E3) ,  and  (E4)  represent  a  simple  configuration 
of  the  beam  which  may  be  useful  in  applications.  Other  combinations 
of  number  density,  axial  and  azimuthal  velocities  are  possible,  how¬ 
ever  (see  (Davidson,  1974;  20,  21)  for  a  general  equation  which  these 
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'.1 

3  b" 
3t 

=  0 

(3.47) 

.i 

A 

j 

r 

3  B  ® 
3t 

=  3 

3r 

(3.48) 

• 

^B® 

3  t  ^ 

=  -1  3  rEg^ 

r  3  r 

(3.49) 

■I 

1 

j 

3  rE® 

3  r 

= 

e 

(3.50) 

3 

Note  that  equation  (3.47)  implies  depends  only  upon  the  value  of 
r  ,  and  that  equation  (3.51)  implies  that  the  form  of  is 


(r  ,  t)  =  - 
r  r 


(r>0) 


where  K  is  an  arbitrary  constant.  IfB^(r,t)is  to  remain  bounded  as 
r-^0  ,  then  K  =  0  and,  thus, 


B^(r,t)  =  0 


(r>0,  t>.0) 


Thus  far  assumptions  (A1)  through  (A6)  have  been  implonented. 
Assumption  (A7)  requires  development  of  a  specific  equilibrium  solu¬ 
tion.  The  "rigid-rotor"  equilibrium  (Davidson,  1974:  30;  Krall  and 
Trivelpiece,  1973!  116,  117)  is  well  known  in  plasma  physics.  A  deri¬ 
vation  of  this  equilibrium  solution  for  a  nonrelativist ic  plasma  is 
now  given.  Use  of  this  specific  equilibrium  is  not  required  in  gen¬ 
eral,  however,  since  any  suitable  equilibrium,  analytically  or  numeri¬ 
cally  derived,  is  suitable  for  the  linearization  process. 
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apply 

assumptions  (A5)  and  (A6).  For  these  reasons 

vector  notation  is 

not  used  here. 

Assumptions  (A5)  and  (A6) ,  uniformity  in  the 

0  and  z  direc- 

tlons, 

respectively, are  now  Invoked  by  neglecting  terms  in  equations 

(3.28)  through  (3.39)  which  involve  3  and  3  .  This  results  in 

38  3  z 

the  following  nonlinear  syston  of  equations  (the  argument  (r,t)  has  been 

dropped  for  notational  convenience): 

3  n 

3  t 

=  -1  3  (nV  ) 
r  3  r 

(3.40) 

3  V 

3t 

2 

=  -V  3  V  +  ^9  +  B®-V  B®] 

V- —  r  —  Fi  r  8  z  z9 

or  r 

+  -^[E^+V  B®-V  B®] 
m  r  9  z  z  8 

(3.41) 

3t  ^ 

=  -V  3  V,  -  ^r'^a  +  ^[E^+V  B®-V  B®] 
r^0  m  ^  zrrz 

dr  r 

+  -^[e^+v  b®-v  b®] 

m  0  z  r  r  z 

(3.42) 

3  V 

3t  ^ 

=  -V  3  V  +  •3-[e®+V  af-VaB®] 

V- —  z  m  z  r0  dr 

0  r 

+  -^[E^+V  B®-Vj  B®1 
mz  r9  9r" 

(3.43) 

--  e" 

T —  r 
t 

=  -  u  „  q  c  “  n 

0  '  r 

(3.44) 

3t  " 

=  -  J  ^  q  c  ’  n  V-j  -  c  '  3  B  ® 

3r  ^ 

(3.45) 

3  t  ^ 

=  -J  qc^nV  +  c'  3  (rB®) 

r  3r 

(3.46) 
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3V  =-V3V  -93V  -V3V 

at  dr  ra:3  dz 


+  ^[E^  +  V  B®-V„B^]  +  ^[E^+V  B®-V^B®] 
m  z  r9  9  r-*  m  z  r9  9  r 


Ur,  V  r  U 

r  9  9  r  • 


3  E  = 
—  r 
3  t 


-uqc^v  +  I^b"-  l-E,"] 

0^  r  r39z  3z^-' 


1_%  = 
o  t 


-u^qc-v^  +  -  Itb!] 


■  o  z  r  o  r  z 


..  2,,  ,  ^  nS  13-S-i 

c  E  =  -'^  qc  V  +  c^L-  — rB^  -  —  J 

ttt  z  0^  z  r3rJ  r39r 

a  t 


3  B  =  0  E,  -  1  3  E' 


3t 


a  z 


39 


=  3  E  -  3  E' 


3r 


o  z 


a_B®  =  1  9_E^  -  I  ^rE® 
3t^  r39'^  rSr 


J.  ^rE®  +  J_  ^E®  +  3  E^ 
r3r  ^  r  99  3z^ 


CL_n 


l3rB‘'  +  l3B3  +3B 
r  0  r  r  00  3  z 


(3.31) 

(3.32) 

(3.33) 

(3.34) 

(3.35) 

(3.36) 

(3.37) 

(3.38) 

(3.39) 


A  far  more  compact  statement  of  these  equations  is  usually  given  in 
plasma  physics  texts  by  the  use  of  vector  notation  (see  (Krall  and 
Trivelpiece,  1973!  85,  86;  Davidson,  1974:  14,  15;  Shkarofsky  et  al, 
1963:  12,  21;  Montgomery  and  Tidman,  1964:  12)).  These  compact  forms, 
however,  can  be  confusing  to  those  not  familiar  with  the  notation. 
Furthermore,  the  detailed  expressions  above  are  needed  in  order  to 
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orthogonal  system  which  has  been  rotated  about  the  z-axis  by  the 
angl e  9 : 


-sin  9 


sin  9 


cos  9 


0 


Letting  [  j  .  ’i,  .  ^  represent  [Vj,V^,V^]'^  , 

'J’  _ 

[  B ^  ,  B  ^  ,  B ^  ]  ,  then,  the  new  functions  [ are  defined  by 


r 

ill 

V0 

= 

X 

—  ~ 

’^2 

N 

1 _ 

1 

u> 

Using  these  definitions,  equations  (3.22)  through  (3.27)  can  be 


expressed  in  terms  of  the  new  unknown  functions  n  ,  v  ,  ,V  ,  E  , 

r  ^  z  r 

s  s  s  s 

E. ,  ,B^  ,  with  independent  variables  (r,9  ,z,t)  (the  argu¬ 

ment  (  £iC)  is  dropped  for  notational  convenience): 


3  n  =  -jAv  )  -  d  nV  -  1  3  nV. 


-  3  nV 


(3.28) 


3  V  = 
—  r 
1 


•V  9  V  -  ^9  3  V  +  -  V  3  V 

or  r  ci9  r  3  z 


+  ^[e'^+V,  B^’-V  B,®]  +  a[E®+V,  B®-V  B  ®  1 
^  r  :z  z9“  -‘•  r  3z  z3^ 

m  m 


(3.29) 


—  9 

r 


^9  3  V  -  -  V  3  V 

-  0:3-  o  -  Z- —  z 

r  o9  r  o  z 


+  a'  El  +V  B'  -V  B  1 
„  -  "  z  r  r  z  ' 


+  q[Ef+V  B®-V  B®] 
j  z  r  r  Z 


(3.30) 
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0 

0 

0 

0 

0 

b“ 


with  explicitly  defined  in  equation  (3.60).  The  operator 

reduces  considerably,  in  this  case,  as  follcws: 


0 

0 

M 

p  0  y 

-b“C3 

.  £ 

-f 

0 

m 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.  — 

The  nonlinear  li  can  be  approximated,  then,  by 
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^  i_i  I  ,  ■  I  iji  wa. 


3.67) 


G (E®, B® ,u(t)  ) 


G(E®’°  .B®’'’  .u®  )  +  ^^0  (3.68) 


Using  the  approximations  of  equations  (3.62),  (3.68)  in  equation 
(3.61)  yields 


^  ii(c)  =  f^(u“  +  iu(t))  ^  KuM  +  G  (E®  ’  “  ,  B®  ’  “  ,u°  ) 


+  11^';  (lii(  t  )  )  +  IG^O  .iB®  ,;^(t  )  ) 


This  implies,  then,  that 


d 

d  t 


iji(  t ) 


il^j0(lu(t))  +  6G^o  (5E^,5B^.5u(t)  ) 


(3.69) 


This  equation  represents  a  linear  approximation  to  the  nonlinear  sys¬ 
tem  of  equations  (3.40)  through  (3.46)  and  (3.48)  through  (3.50).  A 
similar  model  can  be  developed  for  the  region  r>R  .  In  general, 
these  two  models  must  be  solved  simultaneously,  and  their  solutions 
must  satisfy  further  mathematical  constraints  at  r=R  .  Consequently, 
only  small  excursions  of  the  actual  beam  radius  from  the  equilibrium 
solution  radius  are  allowed  before  a  relinearization  must  be  per¬ 
formed. 


Linear  Model.  To  summarize  the  development  thus  far,  the 
expressions  derived  in  equations  (3.64)  and  (3.67),  based  on  the  rigid 


rotor  assumptions,  are  substituted  into  equation  (3.69)  to  yield 


5E^(r,t)  -  cjr6B^(r,t)  -  V^cB^Cr.t) 

6Eg®(r,c)  +  V°5B®(r,t) 

6E^(r,t)  +  ajr£B^(r,t) 

0 
0 
0 
0 
0 

An  appropiate  Initial  condition  for  this  system  of  equations  is 
any  "small"  perturbation  from  the  equilibrium  solution  ^1°  #  Let  ja 
be  any  initial  condition  for  the  system  of  equations  (3»6l)  such  that 

l|u-u°||^  <  d 

where  d  is  a  sufficiently  small  positive  constant.  The  corresponding 
initial  condition  for  equation  (3.71)»  for  any  such  then,  is 


^(t)  = 


m 


I 


M 


given  by 


The  Electrostatlo  Approximation  Model 


An  additional  assumption  simplifies  the  single  degree  of  freedom 
model  CO  nsi  der  ably : 

A8.  perturbed  self  magnetic  field  effects  are  negligible 

This  assumption  is  occasionally  invoked  in  the  study  of  plasmas  and  is 
referred  to  as  the  "electrostatic  approximation"  (Davidson,  1974:  42), 
The  linear  model  implied  by  assumptions  (A1)  through  (A7)  and  the 
electrostatic  approximation,  assumption  (A8),  is  now  developed. 

Assumption  (A8)  eliminates  the  equations  for  WgCt)  and  Wg(t) 
(see  equation  (3.71))  immediately  since  these  represent  (approxima¬ 
tions  to)  the  time  derivatives  of  6:^®(r,t)  and  6B®(r,t),  respec¬ 
tively.  Furthermore,  the  equations  for  Wg(t)  and  W7(t)  are  elimi¬ 
nated  in  the  electrostatic  approximation.  This  is  Justified  as  fol¬ 
lows.  Consider  equation  (3.28)  for  the  perturbed  electric  field 
^E®(x,t)  : 


i  i  k  3  X  . 
J 


5Ef(x,t)  =  (x  ,  t )  =  0 


k 


1  — 


(i  =  1,2,3) 


Since  the  curl  of  the  perturbed  field  vanishes,  it  must  be  expressible 
as  the  gradient  of  a  scalar  field  (Marion,  1965:  105-108): 

cE.(x,t)  =  -|^i(x,t)  (i  =  1,2,3) 

i 

Expression  of  this  result  in  the  cylindrical  coordinate  system  and 
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Assumptions  EL  through  E5  and  A1  through  A8  have  been  introduced 
in  order  to  develop  a  sufficiently  simple  linear  model  for  a  demon¬ 
stration  of  semigroup  theory  techniques.  These  techniques  are  not 
dependent  upon  the  numerous  assumptions  invoked  above.  The  rigid 
rotor  equilibrium  is  merely  one  of  many  equilibrium  solutions;  in 
fact,  numerically  generated  equilibria  can  be  developed  for  a  differ¬ 
ent  set  of  assumptions  than  El  through  E5.  The  fundamental  problems 
to  be  addressed  in  any  infinite-dimensional  system  remain  the  same, 
however.  One  must  demonstrate  that  the  model,  with  appropriately 
chosen  function  spaces,  is  well  posed.  This  is  equivalent  to  showing 
the  operator  A  in  the  abstract  Cauchy  problem  (equations  (2.1), 
(2.2))  is  in  one  of  the  C-classes  described  in  Chapter  II.  In  Chap¬ 
ter  IV,  the  electrostatic  approximation  model  is  used  to  illustrate 
the  theory.  Specifically,  appropriate  spaces  are  chosen  and  an  ana¬ 
lytic  solution  is  obtained.  In  more  complicated  models,  numerical 
methods  will  generally  be  required  a.  d  determination  of  appropriate 
spaces  will  undoubtedly  be  more  difficult,  but  the  basic  principles 
remain  the  same. 


Cone lus ion 

Models  of  the  dynamic  behavior  of  a  charged  particle  beam  have 
been  developed  which  vary  widely  in  complexity.  The  most  accurate 
models  consist  of  the  microscopic  descriptions  and  involve  six  inde¬ 
pendent  space-like  variables.  Macroscopic  descriptions  are  less  com¬ 
plex  and  involve  at  most  three  independent  spatial  variables,  but 
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require  additional  assumptions.  Linearizations  of  both  types  of 
descriptions  yield  models  with  the  abstract  Cauchy  problem  structure. 

A  single  degree  of  freedom  linear  perturbation  model  has  been 
developed  based  on  a  physically  reasonable  equilibrium  solution.  This 
model  is  novel  in  that  it  incorporates  the  effects  of  the  external 
fields  as  controls,  and  it  is  expressed  as  an  abstract  Cauchy  problem. 
As  a  result,  a  new  particle  beam  model  is  now  available  to  researchers 
in  the  control  community  in  a  form  which  is  directly  useful  for  fur¬ 
ther  analysis. 

The  final  result  of  this  chapter  has  been  the  development  of  a 
particularly  simple  model  which  has  a  closed-form  solution.  It  is 
valid  whenever  the  perturbed  self  magnetic  field  effects  are  negligi¬ 
ble,  and  this  is  the  situation  so  long  as  all  perturbed  velocity  com¬ 
ponents  are  sufficiently  small.  A  solution  of  this  model  is  developed 
in  the  following  chapter. 
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IV.  Analysis  of  the  EleGtrostatlc  Approximation  Model 


Introduction 

An  analysis  of  the  differential  equations  in  the  electrostatic 
approximation  model  is  now  presented.  At  the  outset,  conventional 
methods  of  classifying  systems  of  partial  differential  equations  are 
discussed.  The  system  of  equations  (3.72)  does  not  fall  into  any  of 
these  conventional  classifications,  at  least  by  most  authors'  defini¬ 
tions.  In  fact,  no  treatment  of  systems  with  this  particular  struc- 
R  ture  could  be  found,  by  this  author,  in  control  theory  literature. 

Consequently,  fundamental  concepts  must  be  applied  to  the  problem  at 
hand. 

To  this  end,  a  simple  example  of  a  system  which  is  similar  to 
that  of  equation  (3.72)  is  introduced.  This  trivial  example  provides 
insight  as  to  how  one  might  choose  an  appropriate  underlying  Banach 
space  for  these  types  of  systems. 

Both  physical  considerations  and  insight  from  this  example  are 
then  used  in  selecting  a  Banach  space  for  the  electrostatic  approxima¬ 
tion  model.  The  matrix  of  operators  in  equation  (3.72)  is  shown  to  be 
the  generator  of  a  strongly  continuous  group  on  this  space.  The  asso¬ 
ciated  semigroup  of  operators  is  then  constructed,  and  a  closed-form 
solution  of  the  homogeneous  abstract  Cauchy  problem  associated  with 
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equation  (3,72)  is  given.  The  nonhomogeneous  solution  follows  immedi¬ 
ately,  in  ll^t  of  Theoran  2.12,  for  a  broad  class  of  inputs. 


Conventional  System  Classifications 


Consider  the  following  system  of  partial  differential  equations: 


t) 


A(x,  t)^u(x,  t) 


+  b  (u  ;  X  ,  t  ) 


(4.1) 


The  unknown  vector-valued  function  assimies  values  in  ,  A  is  an 
C  h 

n  -order  matrix-valued  function  of  (x,t),  and  ^  is  a  (possibly  non¬ 
linear)  function  of  ^  as  well  as  (x,t).  Systems  of  equations  with 
this  structure  are  sometimes  classified  as  hyperbolic,  parabolic  or 
elliptic. 

Equation  (4.1)  is  called  hyperbolic  at  the  point  (x,t)  if 
(1)  all  roots  of  the  polynomial  P(A;x,t),  defined  by 

P(X;x,t)  =  -  A(x,t)] 

are  real  and  (2)  if  there  exists  a  full  set  of  linearly  independent 
eigenvectors  (Courant  and  Hilbert,  1962:  425;  Garabedian,  1964:  96). 

Some  authors  prefer  to  define  system  (4.1)  as  hyperbolic  only  if  the 
polynomial  P(X;x,t)  has  n  distinct  roots,  while  others  refer  to 
such  a  system  as  "strictly  hyperbolic"  or  "hyperbolic  in  the  narrow 
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sense"  ( Zaclunonoglou  and  Thoe,  I976:  362).  With  this  minor  exception, 
there  is  good  agreement  among  authors  of  partial  differential  equa¬ 
tions  texts  on  the  definition  of  hyperbolic  systems.  If  the  matrix 
A(x,t)  is  symmetric,  then  it  is  well  known  that  a  full  set  of  li&- 
early  independent  eigenvectors  exists.  Consequently,  various  treat¬ 
ments  of  equation  (4.1)  have  been  undertaken  under  this  simplifying 
assumption  (Russell,  1978:  647;  Fattorini,  1983:  146). 

There  is  far  less  agreement  on  the  definition  of  a  parabolic 
system,  however.  Hellwig  (1964:  70)  defines  system  (4.1)  to  be  par¬ 
abolic  if  the  polynomial  P(X;x,t)  has  precisely  v  distinct  real 
roots,  where  l^v^n-1  .  few  authors  aillcw  such  a  broad  definition, 
however.  Various  restrictions  are  usually  imposed  on  equations  with 
the  structure  of  (4.1)  in  order  to  preserve  some  of  the  properties  of 
scailar  parabolic  equations  (e.g. ,  the  heat  equation)  (Eidel'man,  1969: 
3).  As  a  result  of  vairious  authors'  viewpoints,  we  have  systems  of 
equations  which  are  defined  as  "parabolic  in  the  sense  of  Petrovskiy," 
"parabolic  in  the  sense  of  Shilov,"  or  "parabolic  in  the  sense  of 
Shirota"  (Eidel'man,  1969:  444-453).  No  universally  accepted  defini¬ 
tion  of  a  parabolic  system  of  partial  differential  equations  has  yet 
emerged. 

While  there  is  general  agreement  on  the  definition  of  an  elliptic 
system,  many  texts  on  partial  differential  equations  omit  any  dis¬ 
cussion  of  such  systems.  Both  Hellwig  (1964:  70)  and  Zachmanoglou  and 
Thoe  (1976:  362)  define  the  syston  (4.1)  to  be  elliptic  at  the  point 
(x-t)  if  the  polynomial  P(A;x,t)  has  no  real  eigenvalues.  Ellip¬ 
tic  systems  do  not  often  arise  in  initial  value  problems  (Courant  and 


and, 


h  (  X  ,  t  ) 


0  0  0  0  0  0 


1  0  0  0  -V' 


-0!  X 


0  1  0  V  “  0  0 


0  0  1  cox  0 


0  0  0  0  0  0 


5  E  ^  ( X  ,  t ) 
5  ( X  ,  t ) 


5  E  ®  ( X  ,  t ) 
z 


0  B  ^  ( X  ,  t ) 


^  (x , t ) 


6  (x  ,  t ) 


Consequently,  the  polynomial  P(X;x,t)  is  independent  of  (x,t)  , 
and  is  given  by 


P(X,t)  =  det(Xl  -  A^)  =  X’ 


It  is  easily  verified  that  there  are  only  four  linearly  independent 
eigenvectors  associated  with  the  eigenvalue  X  =  0,  Hence,  this  sys- 
tan  is  neither  hyperbolic  nor  elliptic.  Furthermore,  it  is  not  para¬ 
bolic  under  any  of  the  definitions  mentioned  above  except  for  Hell- 
wig's.  Unfortunately,  in  contrast  to  the  hyperbolic  case,  no  exten¬ 
sive  treatments  of  systons  of  this  type  have  been  found  in  control 
literature,  so  equation  (4.2)  will  be  analyzed  from  fundamental  con¬ 
cepts.  To  this  end,  a  simple  example  of  such  a  parabolic  system  is 


useful. 
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An  Illustrative  Example 


Consider  the  two  coupled  partial  differential  equations 
|^Uj(x,t)  =  -|^U2(x,t)  (4.3) 

2(x,t)=0  (4.4) 

with  initial  data 

Uj (x,0)  =  u°  (x)  (4. 5) 

u^ (x,0)  =  u°  (x)  (4. 6) 

for  O^x^l,  t^O .  Equations  (4.3),  (4.4)  have  the  structure 

of  equation  (4.1)  with  ^(£;x,t)  =  0,  and 

A(x,t)  =  A  = 

Note  that  A  has  the  single  eigenvalue  X  =  0,  and  there  is 
but  on  linearly  independent  eigenvector  associated  with  this 
eigenvalue.  The  following  two  propositions  are  now  proven: 

Proposition  1:  The  abstract  Cauchy  problem  of 

equations  (4.3)  -  (4.6)  is  not 
well  posed  for  Xj=  X2  =  L^(0,1). 

Proposition  2:  The  abstract  Cauchy  problem  of 

equations  (4.3)  -  (4.6)  is  well 
posed  for  Xj  =  L^(0,1), 

X,  =  (0, 1). 


d 

dx 


Proof  o f  Proposition  1 
Analysis  of  the  operator 


A  = 


0 

0 


reveals  that,  for  z  0 , 


(zl  -  = 


I 

z 


0 


1  d 
z^  dx 


1 

z 


Now  (zl  -  A)~^  is  not  in  the  set  B(X)  for  any  value  of 
since,  for  the  specific  choice 


X 


one  has  that 


(zl  -  A)  "^^(x)  = 


However ,  x  i  Xj  =  L^(0,1),  and  so  (zl  -  A)“^  ^B(X). 

Consequently,  the  spectrum  of  A  is  the  entire  complex 
plane  and,  thus,  A  ^  G(M,8)  for  any  pair  (M,B)  (recall 
Definition  2.3,  page  II-IO).  By  Theorem  2.2,  then,  the 
abstract  Cauchy  oroblem  (4.3)  -  (4.6)  is  not  well  posed 


of  the  Eleotroetatic 


Model 


The  electrostatic  approximation  model  developed  in  Chapter  III 
(see  equation  (3.72))  is  a  nonhomogene ous  abstract  Cauchy  problem: 


where 


=  A  w  ( t )  +  g  ( t ) 


w(0)  = 


0  a^j 


(4.7) 


(4.8) 


The  symbols  represent  constants  in  the  matrix  of  operators  in 

equation  (3.72),  and  the  operator  is  defined  on  page  III-31.  The 
mapping  w  takes  values  from  the  nonegative  real  line  into  a  function 


space  X  =  n  X  .  —  i.  e. , 


L’V' 


5n (x , t ) 


6V^(x,t) 


w(t)=  6v„(x,t)  eX  (t^O) 


6V^ (x,t) 


5  (x , t ) 


Analysis  of  the  illustrative  example  above,  equations  (4.3)  -  (4.6), 

suggests  that  the  choice  of  X  has  a  profound  effect  on  the  well- 
posedness  of  such  systems.  Specification  of  the  is  now  made  based 
upon  both  the  physics  of  the  problem  and  mathematical  insist  obtained 
from  the  example. 

If  n(r,t)  denotes  the  number  density  in  a  cylindrical  beam  of 
radius  R  with  axial  and  azimuthal  symmetry,  then  the  total  number  of 
particles  in  a  unit  length  of  the  beam  at  time  t,  N(t)  ,  is  given  by 


N(t)  =  27r 


/n(.. 


t)rdr 


The  number  density  can  be  expressed  as  an  equilibrium  value,  n°  (r)  , 
plus  a  perturbed  number  density  6n(r,t)  ,  and  thus. 


N(t)  =  2it 


rR  r  R 

/  n°(r)rdr  +  2tt  J  6n(r,t)r< 


One  natural  choice  for  the  norm  of  the  perturbation  5n(r,t),  then, 
is  the  L  ^  ( 0  ,  R )  norm  of  r5n(r,t); 


l|r5n(r,t)||Li  = 

I  j  5n(r,t)r| 

0 


Define  (0,R)  to  be  the  space  consisting  of  all  functions  g  for 
which 


(0,R) 


(O.R) 


\9 


In  Appendix  C  it  is  proven  that  (O.R)  is  a  Banach  space.  The 

space  X^is  now  defined  to  be  M^(0,R). 

The  Sobolev  space  (0,R)  is  selected  for  the  spaces  Xj 
through  Xj,  Unlike  the  choice  for  Xj  ,  this  selection  is  motivated 
more  by  the  results  of  the  example  in  the  last  section  than  by  physi¬ 
cal  considerations. 

Since  the  spaces  Xj  through  Xj  chosen  above  are  all  Banach 

5 

spaces,  the  Cartesian  product  ^  =  II  X^  becomes  a  Banach  space. 

i=  1 

For  convenience,  define  the  norm  of  the  space  X  by 

I 
\ 
! 
( 

I 

•!  i 


Nllx 


max 

i=  1  ,  .  .  .  5 


Nil 


With  this  choice  for  X  the  linear  operator  A  in  the  electrostatic 


approximation  model  can  be  shown  to  be  a  bounded  operator 


Proof 

In  order  to  show  that  AeB(X)  ,  it  is  sufficient  to  show 
that 

ll^ill  <  Kill II 

for  some  K  >  0  (see  page  II-3).  Denoting  by  i-te- 

qualities  for  ||^|j  eire  derived  as  follows; 


Consequently, 


IlSilIx  i  Ci2l|f,llx. 


Il®2  il  X 


II82IIX2  =  11^23^3  5)1 


By  the  triangle  inequality,  then,  and 
since  X  =X,  =  x,  =  X^, 

2  3  4  5 


X2 


P2  X  =  '^23*  ^3  X,  '^24!  ^^4  X  1^25!  ^5|IX 


The  middle  term  on  the  ri^t  hand  side 
obeys  the  following  inequalities: 


xf.llv  =  /  x^f^(x)^  +  [D^(xf  ^  (x)  )  ]^dx 

4 


■/ 

0 

./: 

0 

r 

<  /  2(x 

^  r\ 


f^(x)2  +  2[fjx)2  +  (xf;(x))']dx 


+  2)[fjx)^  +  f;(x)^]dx 


f 


<  2  (R"+2)  /  f  ^  (x)^  +  f;  (x)Mx 


(R’+2)  l|f 

^  2V 


and  so 


Mix  <  [2(r2  -+  -2)]-^  jlf 


X, 


1 

j 

■’j 
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As  a  result, 


v''2  (R’  +  2) 


+ 


a . 


+ 


2  4 


+  c 


slhllx 


a  X  f 

2  4 


!  J/2(R2+2)  llf^ 


(The  argument  is  identical  to  that  for 
the  term  above.) 


4 


From  these  inequalities,  it  can  be  seen  that  the  norm  of 
satisfies  the  following: 


The  constant  K  is  given  by 


K  = 


max 
i=  1  ,  .  . 


and  the  C^.  are  either  defined  as  above  or  are  zero  if  not 
previously  defined.  Since  A  is  linear  and  ||'^Xiv=^ 

A  £  B(X)  .  f 


Proof  of  the  following  theorem  is  Immediate  in  li^t  of  this 
lemma  and  Theorems  2.2  and  2.6: 


Theorem  4 . 1 

If  the  linear  operator  A  and  the  Banach  space  X  are 
defined  as  above,  then  the  homogeneous  abstract  Cauchy  prob¬ 
lem  associated  with  equation  (4.7) 


=  Aw(t)  (t>0) 


=  w  ■' 


w  ( 0  ) 
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1 


r1 


1 

j 

-j 


kMl 


is  well  posed 


Also  from  Theorem  2.6,  the  solution  of  the  homogeneous  problem  is 
simply 

w(t)  =  S(t)w*’ 

where  S(t)  is  the  strongly  continuous  group  generated  by  A  ; 


S  (t) 


n 


j=0 


t  j  A  j 
j  ! 


A  closed-form  expression  for  S(t)  is  now  developed. 

Computation  of  S(t)  involves  the  development  of  general  terms  in 
the  infinite  series  indicated  above.  The  work  is  simplified  by  the 


use  of  a 

parti  tioned 

matrix 

expression  for  A  ; 

» 

! 

0 

-n  '  D 

:< 

0 

0 

0 

i 

! 

0 

0 

a  ,  3 

a  2  u  X 

^2  5 

1 

C 

F 

0 

A  = 

c 

-  a 

-  ? 

0 

0 

0 

= 

0 

0 

G 

1 

0 

-  a  ^  X 

0 

0 

0 

0 

H 

0 

0 

r  a 

0 

0 

0 

Writing  out  the  first 

few  terms  of  A 

one  obtains 
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noisy  environment. 


Parameter  estimation  should  also  be  pursued  since  some  quantities 
in  the  model  are  not  likely  to  be  known  with  great  certainty  (e.g, 
V  ,  -  ,  and  In  fact,  the  relatively  new  method  known  as  mul¬ 

tiple  model  adaptive  control  (Maybeck,  1982:  253)  could  prove  useful 
as  well. 

Full  Linear  Macroscopic  Model.  An  examir,ation  of  the  linear 
macroscopic  model  of  equations  (3.66)  and  (3.67)  reveals  that,  like 
the  electrostatic  approximation  model,  this  system  of  nine  partial 
differential  equations  is  classified  as  a  parabolic  system  under  Hell- 
wig's  defirdtion.  Writing  equation  (3.66)  in  the  form  of  equation 
(4.1),  one  has 


.\  (,  A  ,  t  )  =  A  = 


and,  thus, 


0  -  n '' 

0  0 

0  0 

0  0 

0  0 

0  0 

0  0 

0  0 

0 


0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 


0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 


0  0 

0  0 

0  0 

0  0 

0  0 

0  0 

0  c 

1  0 


0  0  0  -  1  0  0 


0 

0 

0 

0 

0 

-  c  ^ 

0 

0 

0 


r  -  A ' 


C'  )  ’ 


It  is  easily  verified  that  only  four  linearly  independent  eigenvectors 
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number  density,  velocity  field,  and  radial  electric  field,  under  the 


assumptions  required  for  the  electrostatic  approximation  model.  It  is 
not  yet  clear  whether  an  arbitrary  state  can  be  attained  through  the 
action  of  allowable  controls.  In  the  language  of  control  theory,  one 
would  like  to  establish  whether  the  electrostatic  approximation  model 
is  controllable,  approximatelv  controllable  (Russell,  1978:  643),  or 

nei  ther. 

A  second  control  problem  is  that  of  synthesizing  a  regulator  to 
maintain  the  equilibrium  solution  when  the  model  is  subjected  to 
unknown  (or  unmodelled)  inputs.  Generating  a  stable  configuration  for 
a  plasma  in  the  laboratory  is  frequently  difficult.  A  regulator  based 
on  the  electrostatic  approximation  model  mi^t  improve  the  situation 
considerably. 

Eventually  a  controller  which  would  enable  changing  the  state  of 
the  beam  to  a  new  equilibrium  might  be  sought.  Adaptive  control  would 
be  necessary  if  the  new  equilibrium  were  to  be  far  from  the  original 
equilibrium. 

Observability  Problem.  A  means  of  detecting  the  state  of  the 
system  is  required  in  order  to  design  a  controller.  Analytical  stud¬ 
ies  of  various  sensor  configurations  can  now  be  performed  with  the  aid 
of  the  electrostatic  approximation  model  and  its  solution. 

Deterministic  studies  are  recommended  first,  to  determine  general 
sensor  chairact eristics  required  in  order  for  the  system  to  be  observ¬ 
able  or  distinguishable  (Russell,  I978:  645).  Once  some  general 
requiranents  of  sensors  are  determined,  stochastic  analyses  should  be 
performed  to  determine  how  well  one  can  estimate  the  beam  state  in  a 
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(see  (Courant  and  Hilbert,  1962:172,173)).  These  models  do  not  fit 
into  most  classification  schemes  for  systems  of  first-order  partial 
differential  equations:  however,  they  appear  to  be  physically  signif- 
ioant.  In  fact,  the  electrostatic  approximation  model  has  a  nonstand¬ 
ard  structure,  yet  a  unique  solution  to  this  system  does  exist  as  is 
shown  in  Chapter  IV.  Perhaps  this  result  Indicates  a  need  for  a  bet¬ 
ter  classification  system  than  presently  exists.  (Semigroup  theory, 
as  applied  to  the  abstract  Cauchy  problem,  may  provide  insight  in  this 
direction. ) 

Contributions  to  the  field  of  plasma  physios  consist  of  (1)  a 
solution  of  the  electrostatic  approximation  model  (Chapter  IV),  and, 
(2)  an  introduction  to  (and  a  demonstration  of)  the  application  of 
semigroup  theory  to  collisionless  plasma  dynamics  problems.  The  solu¬ 
tion  of  the  electrostatic  approximation  model  is  a  closed-form  solu¬ 
tion  and  has  not  appeared,  evidently,  previously  in  the  literature. 
It  describes  the  electromechaniccil  oscillations  of  a  very  simple  beam 
dynamics  model.  Under  certain  approximating  conditions,  the  beam  is 
shown  to  oscillate  at  the  plasma  frequency  as  one  might  expect.  The 
full  potential  of  the  techniques  employed  herein  has  only  begun  to  be 
realized  in  this  area  of  plasma  physics. 


Suggested  Areas  for  Further  Research 

Control  Problem.  Equation  (4.10)  provides  an  explicit  means  of 
predicting  the  effects  of  external  electric  and  magnetic  fields  on  the 


V.  Summary  and  Suggested  Areas  for  Further  Research 


Summary  of  Research  Results 

Significant  contributions  have  been  made  in  this  research  effort 
to  three  distinct  fields:  (1)  control  theory,  (2)  applied  mathematics, 
and  (3)  plasma  physics.  These  contributions  are  now  briefly  dis¬ 
cussed. 

The  single  most  significant  accomplishment  in  this  research  is 
the  laying  of  a  foundation  for  application  of  modern  control  theory 
techniques  to  the  beam  dynamics  problem.  This  foundation  consists  of 
three  separate  blocks.  First,  a  concise  description  of  relevant  semi¬ 
group  theory  results  is  given.  Secondly,  a  full  spectrum  of  beam 
dynamics  models  is  developed.  Finally,  a  specific  model  has  been 
exploited  which  fully  illustrates  semigroup  theory  techniques.  The 
closed-form  solution  of  this  model,  with  external  controls  included, 
is  in  itself  significant,  but  more  importantly,  the  solution  process 
used  serves  as  a  pattern  for  future  control  theory  research  efforts. 

Two  aspects  of  this  research  are  of  interest  to  applied 
mathematicians.  Development  of  semigroup  theory  into  a  useful  tool 
requires  more  documented  accounts  of  actual  applications:  this  report 
represents  one  additional  such  account.  Another  significant  result  of 
interest  to  applied  mathematicians  is  the  form  of  some  of  the  systems 
of  PDE  in  Chapter  III,  The  structure  of  some  of  the  models  therein  is 
neither  totally  hyperbolic  nor  hyperbolic  in  the  more  general  sense 
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number  density,  velocity  field,  and  radial  electric  field  evolve  in 
time  from  an  arbitrary  initial  condition,  but  it  predicts  their  evolu¬ 
tion  in  the  presence  of  external  fields  as  well.  In  a  broader  con¬ 
text,  by  using  solution  techniques  that  involve  elements  of  the  semi¬ 
group  theory  of  operators,  this  powerful  and  elegant  theory  is  now 
more  accessible  to  both  plasma  physicists  and  control  theory  research- 


definition,  but  no  extensive  techniques  for  solving  such  systems  exist 
in  the  literatxnre.  By  an  application  of  semigroup  theory,  and  by  a 
careful  selection  of  the  underlying  Banach  space,  however,  this  model 
was  trainsformed  into  a  well-posed  abstract  Cauchy  problem  and  a 
closed-form  solution  was  derived. 

One  significant  conclusion  that  can  be  drawn  from  these  develop¬ 
ments  is  that  systems  of  linear  partial  differential  equations  can  be 
successfully  analyzed  by  semigroup  theory  techniques  regardless  of 
their  conventional  classification.  Various  researchers  in  this  field 
have  recognized  this  fact.  For  example,  Pazy  (  1983;  105,  110)  classi¬ 
fies  equations  of  the  form 

■^w  (t)=Aw(t)  (t^O) 

as  either  hyperbolic  or  parabolic  depending  upon  whether  A  generates 
a  strongly  continuous  semigroup  or  an  analytic  semigroup  (defined  in 
(Pazy,  1983:  60)),  respectively.  Also,  Fattorini  (1983;  173)  classi¬ 
fies  equations  with  this  structure  as  abstract  parabolic  if  every  gen^ 
eralized  solution  of  the  system  is  continuously  differentiable,  and  he 
relates  this  to  the  analytic  nature  of  the  semigroup.  The  electro¬ 
static  approximation  model  provides  a  concrete  example  of  the  need  for 
a  classification  scheme  which  is  based  on  the  properties  of  the  opera¬ 
tor  A  in  relation  to  the  underlying  Banach  space  X  . 

Several  contributions  to  both  plasma  physics  and  control  theory 
have  been  made  in  this  chapter.  In  the  narrower  sense,  development  of 
a  closed-form  solution  of  the  electrostatic  approximation  model  is 


significant  in  Itself, 


This  solution  not  only  describes  how  the 


brackets  is  essentially  unity  and  the  variation  of  C  with  x  can  be 
ignored.  Supposing  this  to  be  the  situation  for  now, 


1 


'2 


In  Brillouin  flew  conditions  (see  the  subsection  "Nonrelativistic 
Rigid-Rotor  Equilibrium"  in  Chapter  III),  ,  and 

u  p 


Q  =  [2cjj^  -  0)^]^  =  oo 
P  P  P 

This  represents  a  limiting  situation  since  the  rigid  rotor  equilibrium 

can  exist  only  ifeo^  >  2oa^  i.e.,  Cl  >  o)  . 

'  c  =  p  ’  p 

On  the  other  hand,  if  for  some  combination  of  x  ,  V°  and  co  , 

z  p 

does  vary  substantially  with  the  spatial  location,  it  can  still  be 
interpreted  as  a  frequency,  but  a  different  frequency  of  oscillation 
would  exist  at  each  point  within  the  beam. 


Conclusion 

Conventional  metho'-’s  of  classifying  systems  of  partial  differen¬ 
tial  equations  were  discussed  early  in  this  chapter.  The  electro¬ 


static  approximation  model  was  shown  to  be  parabolic  under  Hellwig's 


with  initial  condition  I  provided  the  corresponding  components  of 
also  satisfy  Poisson's  equation; 


D  = 


0  =  5_,,o 


s"' 


Consequently,  the  initial  condition  vector  w'’  is  constrained  hy  phys¬ 
ical  CO  ndi  deratio  ns  to  lie  in  the  subspace  Q  defined  by 


Q  =  {w'^eXiD  w°  =  -^  wM 

X  5  1 

The  Freouenev  fi.  In  the  expression  for  the  semigroup  S(t)  , 
equation  (i|.9)»  the  symbol  is  seen  to  appear  frequently.  The  phys¬ 
ical  significance  of  fi  is  now  discussed. 

In  the  derivation  of  S(t)j  the  symbol  ^  was  defined  as 


Using  the  definition  of  and  some  algebra  one  has 


For  any  reasonable  combination  of  x,  V°  and  the  factor  in 

Z  P' 


IV-24 


For  some  applications,  the  requirement  for  £.(t)  to  be  continuously 
differentiable  is  too  strict.  Meaning  can  be  given  to  the  expression 
on  the  rigbt  hand  side  of  equation  (4.10)  for  a  much  broader  class  of 
inputs.  For  example,  even  if  g.  is  in  the  set  (0,T)  ,  this  expres¬ 
sion  is  termed  a  "weak"  or  "mild"  solution  (Pazy,  1983:  108;  Fattor- 
ini,  1983:  89). 

f:nniment<i  on  the  Solution 

Two  additional  topics  concerning  the  solution,  the  effect  of 
Poisson’s  equation  on  the  initial  conditions,  and  the  physical  signif¬ 
icance  of  the  frequency  ^  ,  are  now  discussed. 

Initial  Conditions.  Although  the  solution,  equation  (4.10),  is 
correct  for  any  initial  condition  vector  w°  in  X  ,  there  is  a  physi¬ 
cal  restriction  on  w“,  Poisson’s  equation  (see  equation  (3.50)) 

D  E®  (x,t)  =  ^n(x,t) 

X  r  t, 

has  been  used  in  the  rigid  rotor  equilibrium  derivation  of  n”  (x)  and 
E^(x)  j  but  the  perturbed  quantities  Wi  (t)  =  6n(x,t  ),  and  Wj  ( t ) 

g 

=  6E^(x,t)  must  satisfy  Poisson’s  equation  as  well: 


It  is  easily  verified  that  components  Wj  (t)  and  W5  (t)  obey  Pois¬ 
son’s  equation,  if  w(t)  is  the  solution  of  the  homogeneous  equation 


Similar  results  follow  for  the  remainder  of  the  elements  of 


S(t)  .  The  complete  expression  for  the  strongly  continuous  group  is 
shown  in  the  figure  on  the  following  page.  The  notation  "(•)"  in  the 

a/ 

top  row  is  used  to  indicate  that  operates  on  the  product  of  each 
initial  condition  component  with  the  expression  within  the  brackets. 
For  example,  the  (1,2)  element  of  S(t)  would  operate  on  w°  as  fol¬ 
lows; 


[S(t)]j2W°(x)  =  -n” D 


s  t  0 


The  solution  of  the  homogeneous  problem  is  now  complete. 

The  Nonhomogeneous  Solution.  Recall  from  equation  (3.72)  that 


£(  t )  = 


0 

1 

0 

0 

0 


0 

0 

1 

0 

0 


0 

0 

0 

1 

0 


0 

0 


-v^ 


-OJX 


cox 


0 

0 

0 


0 

0 

0 


6  E  (x  ,  t ) 
r 


SE^iXyt) 


6E  (x,t) 
z 


6  (x  ,  t ) 
r 


5  Bg  (  X  ,  t  ) 


6  b  ( X , t ) 

z 


By  Theorem  2.12,  if  the  external  fields  are  such  that  ^(t)  is  contin¬ 
uously  differentiable  for  all  t£[0,T],  then  the  nonhomogeneous  equa¬ 
tion  (4.7)  has  the  solution 


w(t)  =  S(t)w'’  + 


/  S(t- 


s)£^(s)ds  (O^t^T) 


(4.10) 


-21 


2k 


S (t)=I+ 


(2k)  ! 


0  0 
0  (GH)^ 


k=  1 


0  0 


FG(HG) 

0 

(HG)*" 


2k+l 


(2k  +  l) 


k=0 


0  F(GH) 
0  0 
0  H(GH) 


The  (1,2)  element  of  S(t)  can  be  summed  as  follows: 


[S (t) 


00 

■E 


k=0 


•E 


k=0 


2k+l 


(2k+l)  ! 


F(GH) 


(2k+l  )  !  ^”^2  3  "  ^^2  4^^  ^25^52^ 


An  examination  of  GH  reveals  that  GH<0  : 


GH  =  -a^3  -  (a^.x)"  +  a,, a 


2  5  5  2 


=  -(OJ  -20)) 
c 


0)^  V° 
P  z 

2c2 


2  2 
X  -  0) 


Letting  -GH ,  then,  one  has 


(GH)^  =  (-1)^ 


With  this  form  for  (GH)  , 


[S  (t)  ] 


-  ^,0-5  V  (i)^^  t 

12  ~  ^  ^  ( 


2k  .  2k+l 


2k+l 


(2k+l)  ! 


k  =  0 


07, 


s  infi  t 


S (t)=I+t 

0 

0 

G 

rr 

0 

GH 

0 

_0 

H 

0_ 

0 

0 

HG 

+ 


3  ! 


0  F(GH) 
0  0 


0 

G(HG) 


4  ! 


0  0 
0  (GH)^ 


FGHG 

0 


_o 

H(GH) 

0  J 

Lo 

0  (HG)fJ 

0 

F (GH) ^ 

0 

0 

0  FG(HG)^ 

5-!' 

0 

0 

G(HG)^ 

6! 

0 

(GH)^  0 

0 

H(GH) ^ 

0 

0 

0  (HG)^  __ 

(Note  that  since  F  is  an  operator,  and  G,  H  are  matrices,  the  order 
of  these  factors  in  the  expressions  is  crucial.)  A  general  term  of 


this  expansion  can  be  seen  to  be 

0  F  (  G  H ) 


^2 


0 


0 


H(GH) 


G(HG) 


3  ~ 


-i  ! 


FG (HG) 


-  1 


(GH) 


2. 

2 


(HG) 


j  =  2,4,6, 


After  some  manipulation,  then,  S(t)  can  be  expressed  as 


are  associated  with  the  eigenvalue  X  =  0  .  The  space  X  must  be 
chosen  with  this  in  mind  if  the  model  of  equations  (3.66)  and  (3.67) 
is  to  be  a  well  posed  abstract  Cauchy  problem.  A  solution  of  this 
model  would  be  beneficial  since,  unlike  the  electrostatic  approxima¬ 
tion  model,  electromagnetic  effects  are  included.  By  comparing  the 
solutions  of  the  two  models,  then,  an  assessment  of  the  shortcomings 
of  the  electrostatic  approximation  model  can  be  made. 

Two  Degree  of  Freedom  Model.  Although  much  can  be  learned  from  a 
single  degree  of  freedom  model,  many  of  the  current  particle  beam  con¬ 
trol  elements  require  at  least  a  two  degree  of  freedom  model.  Analy¬ 
sis  of  the  dynamic  behavior  of  a  beam  inside  a  quadrupole  with  vari¬ 
able  magnetic  field,  for  example,  could  suggest  totally  new  means  of 
beam  control.  It  is  recommended  that  an  equilibrium  solution  be 
sought  for  a  beam  with  assumption  (A6),  the  axial  symmetry  assumption, 
removed. 

System  Classification.  As  was  mentioned  in  Chapter  IV,  con¬ 
ventional  schemes  for  classifying  systems  of  partial  differential 
equations  are  not  well  suited  to  control  applications.  It  is  recom¬ 
mended  that  further  investigation  be  conducted  to  establish  classifi¬ 
cations  of  such  systems  based  upon  both  the  operator  A  in  the 
abstract  Cauchy  problem,  equation  (2.1),  and  the  underlying  Banach 


space  X 


Symbols  that  are  frequently  used  in  Chapter  II  are  summarized 


below.  Page  numbers  are  given,  following  each  definition,  to  indicate 
where  the  symbol  was  first  introduced  in  the  text. 


f  :  A-^B 


V(f),  R(f) 
R,  C 


Re  (z ) 


I ,  Q 

n 

n  X. 
i=i " 

lP(Q),  H'i(a) 


X 


(u,  V) 


a  function  f  which  maps  each  element  of  set  A 
into  an  element  of  set  B  (II-1) 

domain,  range  of  the  function  f  (II-1) 

real,  complex  field  of  scalars  (II-2) 

real  part  of  complex  number  z  (II-IO) 

re-1,  complex,  n-fold  Cartesian  product  of  R ,  C 
(II-2) 

intervals  in  R,  *^"(II-2) 

Cartesian  product  of  spaces  X.,  i  =  1,2, ...n 
(II-2)  ^ 

Lebesgue,  Sobolev  space  of  order  p  ,  q  over 
interval  ^  (II-2) 

norm  of  f  on  space  X  (II-2) 

inner  product  of  u ,  v  (11-24) 


B(X,Y),  B(X) 


set  of  all  bounded  linear  operators  from  X  into 
Y  ,  from  X  into  X  (II-3) 


Qj|a3 


C(X,Y),  C(X)  set  of  all  closed  linear  operators  from  subset  of 

X  into  Y  ,  from  subset  of  X  into  X  (11-5,6) 

CO 

Cq  (fi)  set  of  all  functions  which  are  continuous,  have 

continuous  derivatives  of  all  orders,  and  which 
have  support  bounded  and  contained  in  Q  (II-6) 


C  [  a  ,  b  ] 


set  of  all  functions  which  are  continuous  in  the 
sup  norm  on  [a, b]  (11-27) 


1  o  c 


set  of  functions  in  (K)  for  every  bounded, 
Lebesgue  measurable  set  K  with  closure  contained 
in  Q  (II-6) 


ordinary,  partial  differentiation  symbols 
(generalized  derivatives  implied  unless  otherwise 
stated)  (11-6,7) 


.  x^n  ,  Where  I  k|  =  Z  k  , 

k=(ki  ,  .  .  •  k  ),  with  k.  nonegative  integers  for 
i=l,...n  (II-7)  ^ 


Gateaux,  Frechet  derivative  of  operator  F  at  x 
(II- 8) 


P  (A) 
a  (A) 

R(z  ,A) 

G  ( M  , ,  C  '  ( M  ,  £  ) 
"support  of  f" 


resolvent  set  of  linear  operator  A  (II-9) 
spectrum  of  linear  operator  A  (II-9) 
(zI-A)~^  where  zep(A)  (Ii-10) 
see  DefirJLtion  2.3  (11-10) 
the  set  of  points;  {x:  f(x)  /  O) 


The  following  symbols  are  introduced  in  Chapter  III  and  are 
summarized  here  for  convenience.  Page  numbers  in  parentheses 
following  the  definitions  indicate  where  the  symbol  was  first  used 
defined. 

m  particle  rest  mass  (III-4) 

q  particle  charge  (III-4) 

c  vacuum  speed  of  li^t  (III-4) 

£(j  absolute  dielectric  constant  (III-8) 

pQ  absolute  magnetic  permeability  of  free  space 

(III-8) 

cyclotron  frequency  (III-27) 
plasma  frequency  (III-27) 

(i,j,k)  =  (1,2,3),  (2,3,1),  or  (3,1,2) 
(i,j,k)  =  (1,3,2),  (2,1,3),  or  (3,1,2) 
otherwise 

iu|  magnitude  of  vector  (III-4) 


microscopic  velocity  vector  (III-5) 


£ 


microscopic  manentum  vector  (III-5) 


n (x , t ) 
V(x,  t) 

P(x.t) 

f 

J (x ,  t  ) 
E (x , t ) 
l(x , t ) 
P 


number  density  (III-5) 

macroscopic  velocity  vector  (III-5) 

macroscopic  momentum  vector  (III-5) 

distribution  function  (III-5) 

/ 

f (x t ) /n (x , t )  n(x,t)  >  0 
i  (Page  III-7) 

0  n  (x  ,  t )  =  0 

current  density  vector  (III-6) 

electric  field  vector  (III-4) 

magnetic  field  vector  (III-4) 

pressure  tensor  (III-6) 

transformation  matrix  (III-19) 


In  Chapter  IV  the  space  (0,R)is  defined  as  follows; 


M>(0.R)  -  <  -) 


^  (0,R)  IKUl^  (0,R) 


It  is  asserted  in  that  chapter  that  (0,R)  is  a  Banach  space  and 
this  is  now  proven. 


rheorem  C.  1 


The  space  (0,R)is  a  Banach  space. 


Proof 


It  is  obvious  that  (0,R)  is  a  normed  linear  space.  To  show 

-  00  1 

completeness,  let  Ig^ljCM  (0,R)  be  a  Cauchy  sequence.  Defining 
functions  f  ^  by 


f^(x)  =  xg^(x) 


(0<x<R) 


it  is  clear  that  (f^}jCL^(0,R)  , 
L  ^  ( 0  ,  R )  si  nee 


But  this  sequence  is  Cauchy  in 


II  m  nil  (0  ,  R) 


^n^llL^  (O.R) 


I|Sm  ShIIm^  (0  ,  R) 

oo 

and,  as  is  Cauchy,  there  exists  an  N  for  every  £>0 


®m”®n  1m^  (0  ,R)  ^ 


Vm , n  >  N 


T  '■ 


(0,R) 


Vm , n  >  N 


1  r  'i  “ 

Since  L  (0,R)  is  complete,  the  sequence  converges  to  an 

element  f  in  (0,R).  Let  a  function  g  be  defined  by 


g (x)  =  -f (x) 


(0<x<R) 


New  g  e  M  ^  ( 0  ,  R )  si  nee 


-  II*  (in  II, 


®‘*M^(0,R)  (0,R)  l•^llL^(0,R) 


OO 

The  sequence  then,  converges  to  g  since 


^IIm^CO.R)  ®HIl^(0,R)  IIlVo  ,R) 


and  f  f  in  the  (0,R)  norm.  Thus,  every  Cauchy  sequence  in 
(0,R)  converges  to  an  element  in  (0,R)  . 


Appendix  12^.  A  Three  Degree-of-Freedom 
Line.ar  Hfi.d.el 


As  stated  in  the  footnote  on  page  III-17.  the  linearization  proc¬ 
ess  is  in  no  way  limited  to  a  single  degree-of-freedom  model.  By 
applying  the  same  techniques  used  in  Chapter  III  to  the  system  of 
equations  on  pages  III-20  and  III-21*  with  the  rigid  rotor  equilibrium 
solution  (equation  3.60),  one  can  derive  a  linear  model  with  the  fol¬ 
lowing  form: 


^w(t)  =  Aw(t)  +  £(t) 


The  symbol  w(t)  is  given  by 


w  ( t )  = 


iSn  (r ,  t) 

6V^(r,t) 

'^•’^0  t ) 

SV^(r,t) 

5E^(r,t) 

6Eg (r , t) 

5  (r  ,  t ) 

5B^(r,t) 

5  Bg  (£,  t  ) 

5  B^  (_r  ,  t ) 


;  r  =  (r  ,  0  ,  z) 


The  operator  A  is  expressed  in  detail  on  the  following  page,  and 


D-1 


^(t)  =  ■^G(r)u(t:) 
^  m  — 


■ 

— 

f  —  1 

0 

0 

0 

0 

0 

0 

(r,t) 

1 

0 

0 

0 

-v° 

z 

-0)  r 

6  Eg  (  I^,  t  ) 

0 

1 

0 

v“ 

z 

0 

0 

■^E®(r  ,t) 

0 

0 

1 

Oir 

0 

0 

5  B  ®  ( £ ,  t ) 

.  .  . 

d  Bg (r  ,  t ) 

0 

5  (  r  ,  t ) 

2  “ 

Note  that  the  form  of  the  differential  equation  above  is  identi¬ 
cal  to  that  of  equation  (3.71)t  the  single  degree-of-freedom  case. 


re 


10 

The  vector  w(t)  is  an  element  of  a  function  space  X  =  ,1  X  whe 

i  =  1  ^ 

each  is  an  appropriately  defined  space  of  functions  defined  on  a 

subset  of  The  dimension  of  the  space  is  now  ten  rather  than 

nine,  since  the  perturbed  radial  magnetic  field  is  identically  zero  in 
the  single  degree-of-f r eedora  case. 
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Dynamic  models  of  a  charged  particle  beam  subject  to  external  electromag¬ 
netic  fields  are  cast  into  the  abstract  Cauchy  problem  form.  Various  appli¬ 
cations  of  intense  charged  particle  beams,  i.  e, ,  beams  whose  self  electromag¬ 
netic  fields  are  significant,  might  require,  or  be  enhanced  by,  the  use  of 
dyramlc  control  constructed  from  suitably  processed  measurements  of  the  state  of 
the  beam.  This  research  provides  a  mathematical  foundation  for  future  engineer¬ 
ing  development  of  estimation  and  control  designs  for  such  beams. 

Beginning  with  the  Vlasov  equation,  successively  simpler  models  of  intense 
beams  are  presented,  along  with  their  corresponding  assumptions.  Expression  of 
a  model  in  abstract  Cauchy  problem  form  is  useful  in  determining  whether  the 
model  is  well  posed.  Solutions  of  well-posed  problems  can  be  expressed  in  terms 
of  a  one-paraneter  semigroup  of  linear  operators.  (The  state  transition  matrix 
for  a  system  of  linear,  ordinary,  first-order,  constant  coefficient  differential 
equations  is  a  special  case  of  such  a  semigroup.)  The  semigroup  point  of  view 
allows  the  application  of  the  rapidly  maturing  modern  control  theory  of  infin- 
ite-dimensioiml  systems. 

An  appropriate  underlying  Banach  space  is  identified  for  a  simple,  but  non¬ 
trivial,  single  degree  of  freedom  model  (the  "electrostatic  approximation  mod¬ 
el"),  and  the  associated  one-parameter  semigroup  of  linear  operators  is  charac- 
.erized. 
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